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1. Prediction by Numerical Methods of Self-Sustained Oscillations in 
a Two-Axis Model of the Nonlinear Continuous-Data LST System 

1-1 . Introduction 

Jt has been demonstrated in [1] that the methods of continuous 
and discrete describing function analysis can be applied to predict 
the existence of self-sustained oscillations in the single-axis model 
of the LST system with nonlinear CMG friction characteristics. 

Furthermore j ^t has been shown in [2] that the stability equations 
as a result^of the describing function analysis may be solved by a 
numerical-iterative technique instead of the usual graphical methods. 

With an appropriate guess of the initial condition, the numerical 

method is found to be quite effective in leading to a convergent solution 

rapidly. 

In the two-axis model of the LST system, the system contains 
two nonlinearities, and the general form of the stability equation in 
the continuous-data case is 

1 + Gg( jio)N(A) + G c (ju)N 2 (A) = 0 (1-1) 

where Gg(ju)) and G^(ju)) denote transfer functions which depend on the 
linear elements of the coupled systems, and N (A) represents the describing 

A A 

function of the CMG nonl inearity. In general, Gg(jw) and G c (jw) are 
functions of frequency w, and N(A) is a function of the amplitude A of 
the assumed sinusoidal input to the nonl inearity. 

The usual graphical method cannot be used to solve for the values 
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of w and A for self-sustained oscillations in Eq. (1-1 ), due to the 
N (A) term in the equation. However, if the system parameters are 

/s 

^ such that the linear term Gg(joo)N(A) dominates over the quadratic 

term, or vice versa, then by neglecting the smaller quantity in Eq. (1-1), 
an approximate solution may still be obtained graphically. 

i\ direct but more time consuming approach of solving Eq. (1-1) would 
be to calculate the terms in Eq. (1-1) for a wide range of values of w 
and A until all combinations which satisfy the equation are found. 

Still 'a third alternative of solving Eq. (1-1) is to use the 

\ 

numerical-iterative method reported in [2]. If an appropriate initial 
guess can b£ made, the numerical-iterative method should lead to the 
exact solutions in an efficient manner. 

All three of the above-mentioned methods have been applied to the 
two-axis continuous-data LST system. It is found that the graphical 
method with approximation and the direct method can provide useful 
information to the final solutions, and the numerical -iterative method 
is effective in arriving at the exact solution on a digital computer. 

The results of these studies are reported in the ensueing sections. 

1 -2. The Continuous-Data Two-Axis LST System Model and Its Stability 

Equation 

Figure 1-1 Shows the signal flow graph representation of the 
continuous-data single-axis LST system. The nonlinear friction of the 
CMG is modeled by the branch with the gain N. If two such models are 
coupled together at the output stage through output torsional coupling, 
the two-axis model of Fig. 1-2 results. This representation of a 
continuous-data two axis LST system may not be a rigorous one from the 
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Figure 1-1. The simplified single-axis continuous-data 
LST system model . 
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Figure 1-2. The simplified two-axis continuous-data LST system model 





5 


structural standpoint. However, the purpose of the present study is to 
develop analytical techniques which are applicable to mul ti-coupl ed 
^ nonlinear systems. It is conjectured that for a two-axis nonlinear 
system, the describing function analysis will generally lead to a 
stability equation of the form of Eq. (1-1). Therefore, for these 
above-lfientioned reasons, the model shown in Fig. 1-2 is considered to 
be adequate. The techniques developed in this report can be applied 
to the prediction of self-sustained oscillations in any continuous- 
data two-axis nonlinear system that is amenable to the describing 
function method. 

With Reference to Fig. 1-2, it can be seen that there are seven 
individual loops, with the following loop gains: 


A, = - — i 

0-2) 

1 s 


NG~ 

A. = - J 

0-3 

z s 

B 1 ■ G 2 G 3 

(1-4) 

b 2 = - g 2 g 3 

0-5) 

C 1 = ' G 1 G 2 G 3 G 4 H 

0-6) 

- G-^GgG^H 

0-7) 

2 2 
D = k‘g‘ 

0-8) 


where G-j , G^, G^, and G^ are transfer functions as defined in Fig. 1-2. 
All system parameters and variables, including the nonlinear describing 
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function N, are consistent with those defined in the single-axis LST 
model reported previously [1,2]. The describing function N(A) is a 
^ function of the amplitude of the input sinusoid of the CMG nonlinearities. 
In the present case, since the two axes are identical, and the couplings 
are symmetric, it is assumed that the amplitudes of the input signals 
at the two nonlinearities are identical. 

The characteristic equation of the coupled system in Fig. 1-2 is 


A = 0. 




where 


A = 1 - (A-j + A 2 + B 1 + B 2 + + C 2 + D) 

+ (A 1 A 2 + A 1 B 2 + A ] C 2 + A-jD + + B-jD 

+ C^ 2 + C.jB 2 + + A 2 D + B 2 D) 

- (A 1 AgD + A 1 B 2 D + B + B 1 B 2 D) 


(1-9) 


0 - 10 ) 


The last equation is simplified if we define 



(1-11) 


Then, Eq. (1-10) becomes 

A = 1 - 2{A + B + C) - D + (A 2 + B 2 + C 2 + 2AB + 2AC + 2BC 
+ 2AD + 2BD) - (A + B) 2 D 


( 1 - 12 ) 


or 
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A = 1 - 2(A + B + C) - D + (A + B + C) 2 - (A + B) 2 D 


Substituting Eqs. (1-2) through (1-8) in Eq. (1-13) yields 


2G«N « « 

A = 1 + ~y~ + 2G 2 G 3 + ZG^O^^H - 


+ — ^ N 2 + G^G^ + (GiG 2 G 3 G 4 H ) 2 + 


2NG 2 G 3 


2NG.G,G 2 G a H , - 2NG,G 2 K 2 

V s — + 2G 1 G 2 G 3 G 4 H - — 

, . .2.,2 n2g 3 G 4 K 2 2NG 2 G 3 G 4 K 2 

2G 2 G 3 G 4 K 2 2 5 


2 2 2 2 
- G 2 G 3 G 4 K 2 


The characteristic equation in Eq. (1-9) can be written as 


A = 1 + G a + GgN + G c rf = 0 


where 


G. = 2G„G- + 2 G.G 2 G.G-H - K 2 G 2 + G 2 G 2 + (G 1 G 2 G 3 G 4 H) 2 

22 22 2222 
+ EG^G^H - 2G 2 G 3 G 4 K 2 - G^G^ 


2G 3 2G 2 G 2 2 Gi G 2 G 2 G 4 H 

b n = + + 

B s s s 

2 2 2 
or 


? ? 
2G 3 G 4 K 2 


(1-13) 


0-14) 


0-15) 


0-16) 


s 


0-17) 
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G„ = -4 - 


r Z r 2 v 2 

G 3 G 4 K 2 


0-18) 


For stability analysis, the characteristic equation of Eq. (1-15) 
is written as 


* 2 

1 + G g N + G c N^ = 0 


where 


A 



'£b 


1 + G, 


A 




1 + G 


A 


0-19) 


0 - 20 ) 


0 - 21 ) 


Equation (1-19) is a function of the system frequency w and the. 
input amplitude A. Since this equation is defined in the complex 
plane, it represents a set of two nonlinear equations with two unknowns. 

When these nonlinear equations have a solution for u> and A, it represents 
a condition of self-sustained oscillations for the system. The oscillations 
may be stable or unstable; therefore, the solutions must always be 
checked for stability. 

1-3. Prediction of Self-Sustained Oscillations By The Approximation 
Method 

It was mentioned in Sec. 1-1 that the graphical method of predicting 
self-sustained oscillations can still be applied to the two-axis LST 
system if one of the last two terms in Eq. (1-1) can be neglected. In 
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other words, the following two conditions may exist: 

1. |G B (ju)N(A)| » |G c (ju))N 2 (A)| 

Then Eq. (1-1) may be approximated by 

1 + Gg( ju)N(A) = 0 0 - 22) 

V 

and the condition of self-sustained oscillations is found from 

the following equation: 

V ju) = ' nW (1 ' 23) 

2. |G B (ju)N(A)| « |G c (ju))N 2 (A)| 

Then Eq. (1-1) is approximated by 

1 + G c (jto)N 2 (A) = 0 0-24) 

and the condition of self-sustained oscillations is found from 

G r (ja>) = - -J— 0-25) 

C N 2 (A) 

Therefore, the graphical solutions involve the plotting of the 
curves for G g ( jou) and - 1/N(A), or G c (jw) and -1/N (A), as the case 
may be. 

Figure 1-3 illustrates the plots of |Gg(jo))N(A) | and |G^( ju))N^(A) | 
versus a> for various values of A, and |N(A)| versus A. The following 
parameters are used for the LST system: 

J„ = 10 5 , = 2.1 , K = 216 , K, = 9700 

V u p 1 



DECIBELS 
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Figure 1-3. 
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H = 600 , K q = 5758.35 , ^ = 1371.02 , K 2 = 100 , K 3 = 3 

Y = 1.38 x 10 7 

The last two parameters, K 2 and K, are the coefficients of coupling 
between the two axes. 

Jhe curves in Fig. 1-3 give information on the ranges of u> and 
A in which the approximation of Eq. (1-22) or Eq. (1-24) is valid. 

Figure 1-4 illustrates the regions in the to versus A plane in which the 
two approximations are valid. The cross-hatched area represents the 
region in which no approximation can be made, and the graphical method 
cannot be used. The criterion of >.20 db and £-20 db is used for 
magnitude comparison for significance. 

The results of Fig. 1-4 show that the graphical approximation 
method is valid for the following ranges of w and A: 

1 . Use 1 + G C N 2 = 0 : 

A £ 10“ 7 to £ 40 rad/sec 

and 60 £ to £ 80 rad/sec 

2. . Use 1 + G„N = 0 : 

A ^ 6 x 10” 5 to £ 600 rad/sec. 

For the region of validity of case 1 above, Fig. 1-5 shows the 
plots' of G c and -1/N 2 . The heavy portions of the curves indicate the 
parts which are valid for stability analysis. Similarly, Fig. 1-6 

A 

shows the case when the equation 1 + G g N = 0 may be used for approxi- 

A 

mation. Again, the intersection between the heavy portions of the G g 
and the -1/N curves would indicate the possibility of self-sustained 
oscillation in the two-axis LST system. Since there are no intersections 
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between the valid curves in Figs. 1-5 and 1-6, the system under study 
is stable for the parameter values used. This result will be sub- 
^ stantiated by the two other methods which are discussed in the following 
sections. 

1-4. A Direct Method of Predicting Self-Sustained Oscillations in 
The Two-Axis LST System 

A direct but time-consuming method of predicting self-sustained 

oscillations in the two-axis continuous-data LST system is to calculate 

all the terms in the following stability equation for a wide range of 
% 

values of w and A: 

v 

V 

1 + Gg( jio)N(A) + G c (jw)N 2 (A) = 0 (1-26) 

Although the direct method would be quite tedious, and it is 
possible that with a given selected increment of variation for the 
values of u> and A, the solution runs may still miss the exact solutions, 
the results in general will give insight to the appropriate guess on 
the initial values for the numerical -iteration method. Therefore, it 
is enlightening to investigate this direct approach at this point. 

/v A O 

Figures 1-7 through 1-9 show the plots of Gg(jw)N(A) + Gq(jw)N (A) 

7 

for y = 1.38 x 10 , and all the system parameters as given previously, 
and for three different sets of values for the coupling coefficients 
l <2 and Kg. Note that for the three cases illustrated, all the trajectories 
for various combinations of A and u) do not intersect the critical point 
which is at 0 db and -180 degrees. Thus the system is stable. However, 
the results show that as the value of is decreased, the trajectories 




Figure 1-5. 




Figure 1-6. 









Figure 1-8. 
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get closer to the critical point. Figure 1-10 shows the trajectories 

when l<2 = 0.1 and = 0.1. It is noticed that several trajectories 

are very close to the critical point. Figure 1-11 shows a magnified 

version of Fig. 1-10 around the critical point. The figure shows that 

-fi -fi -fi 

the trajectories for A = 5 * 10 , 10 , 5 * 10” , are all very close to 

the critical point. This gives indications that there may be more than 
one solution. It will be shown in the next section by the numerical - 
iterative method that there are indeed two solutions for the stability 
equation. One is at A = 5.9867 * 10”^ and u> = 1.88 rad/sec, and the 

i _c 

other is at A - 5.07397 x 10 and w = 4.1086 rad/sec. These solutions 
must still 6e checked for stable or unstable equilibrium solutions. 

1-5. Exact Solution of the Stability Equation by Numerical -Iterative 
Techniques 

In this section the stability equation developed in Sec. 1-2 is 
solved numerically for its exact solutions. The numerical method 
utilized has been described in [2] and is found to be quite effective 
for the two-axis LST system. 

The stability equation, Eq. (1-19), can be written as 

1 + Gg( ju)N(A) + G c (j M )N 2 (A) = 0 (1-27) 


Define 


Gg(joj) - Gg-| + jGji 


(1-28) 


Gj.(ju) = Ggg + jGjg (1-29) 


N(A) = N r -| + jNji 


(1-30) 




Figure 1-10. 
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N 2 (A) = N R2 + jN I2 (1-31) 

where G^-j > 5 ^12* ^Rl * ^1] * ^R2* ^12 sll real 

quantities. 

If Eqs. (1-28) through (1-31) are substituted into Eq. (1-27), 
it becomes 

1 (Gf^-j * ^Il^^Rl = ^ (1-32) 

When the real and imaginary parts are separated, Eq. (1-32) yields 

\ 

A = A r + jAj = 0 (1-33) 

with 


A R 1 + G R1 N R1 ' G I1 N I1 + G R2 N R2 " G I2 N I2 ~ 0 0-34) 

and 

A I = G I1 N R1 + G R1 N I1 + G I2 N R2 + G R2 N I2 = 0 (1-35) 

Further simplification is possible in Eqs. (1-34) and (1-35) 
if we recognize that N R 2 + j N j ^ = (N^ + jN^) but is not necessary 
since the solution of Eqs. (1-34) and (1-35) will be performed on a 
digital computer. 

Equations (1-33) and (1-34) represent two equations in the two 
variables u> and A. As before [2], let us define 

r « i 


x 


A 


) 


(1-36) 
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Then, Eqs. (1-34) and (1-35) can be written as 

F(x) = 0 (1-38) 

The Newton-type quadratically convergent numerical method 
described in [2] can now be directly applied to this two-variable 
system. * 

'ir 

To initiate the iterations, an initial solution is needed. The 
results of the direct calculation method, presented in the previous 
section provide an adequate initial solution. As in that section, the 
parameters which need to be varied are « 2 and K^. Once a solution for 
a given value of « 2 and K 3 is known, and/or l< 3 could be changed 
slightly and the new solution is obtained by using the old solution 
as the initial guess for the new system. 

With l<2 and K 3 set to very small values, the system is decoupled, 
and the solution of the single-axis LST system could be used as an 
initial guess for that of the two-axis LST system. 

If 1< 2 = 0.1 and K 3 = 0.1, the results from the direct calculation 
method shown in Figs. 1-10 and 1-11 indicate that two possible solutions 
exist. The approximate amplitudes and frequencies of the solutions are 

Solution No. 1 o> - 5 rad/sec 

A = 5 x 10' 6 rad 

Solution No. 2 u> = 2 rad/sec 

A = 1 x lCf 6 rad 
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With these approximate solutions as initial conditions, the results 
of the numerical iteration procedure are obtained and tabulated in the 
first set of iterations in Fig. 1-12 and Fig. 1-13. These iterations 
indicate that the actual solutions are 

Solution No. 1 w = 4.1086 rad/sec 

A = 5.0739 x 10‘ 6 rad 

Solution No. 2 a) = 1.88 rad/sec 

A = 5.99 x 10~ 7 rad 

Figures 1-12 and 1-13 also show the change in the solutions when 
K 2 and K 3 are" reduced to K 2 = K 3 = 0.01 and « 2 = K 3 = 0.001, respectively. 
As mentioned previously, the initial guess for a new set of values of 
K 2 and K 3 is the solution with the previous values of « 2 and K 3 . 

With K 3 fixed at 0.001, if K 2 is increased (Figs. 1-14 and 1-15), 
the solutions will change, until K 2 = 8 for solution no. 1 and K 2 =7 
for solution no. 2, beyond which no solution is obtained. It is 
interesting to note that the two solutions move closer together with 
increasing K 2 and almost merge into each other before disappearing 
altogether. The result is analogous to lowering the G(ju>) curve or 
raising the -1/N curve in the single-axis case. 

It can be concluded that for K 3 = 0.001, K 2 should be less than 
8 for a sustained oscillation to occur. In fact, with f< 2 = 8 (solution 
1) and K 2 = 7 (solution 2) the roots of the stability equations do not 
appear to have converged adequately, and well defined roots occur only 
for lower values of A maximum of 20 iterations are attempted before 

the numerical process is terminated. 

If « 3 is kept fixed at 0.2 and IC, is varied, the iterations in 
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Figs. 1-16 and 1-17 show that Kg can be increased only up to Kg = 3 
for solutions to be obtained. 

On the other hand, Kg can be kept fixed and Kg varied; these results 
are shown in Figs. 1-18 through 1-25. In Figs. 1-18 and 1-19 Kg is 
fixed at 0.001, and the maximum value of Kg for a solution to the stability 
equation is approximately Kg = 0.4 for solution no. 1 and Kg = 0.3 for 
solution no. 2. With Kg = 0.1, 1.0 and 5.0, the maximum values of Kg are 
0.4 (solution 1 and 2), 0.3 (solution 1 and 2) and <0.1 (solution 1 and 2), 
respectively. These results are shown in Figs. 1-20 through 1-25. 

Table 1-1 summarizes all the results. In the table, the solutions 
for Kg = Kg t = 0 are the solutions with the single-axis LST system. It is 
seen that when Kg and Kg are very small, the solutions of the two-axis 
system approach those of the single-axis system. 

The data in Table 1-1 can be used to plot the region in the Kg - Kg 
plane in which sustained oscillation can exist. Figure 1-26 shows this 
region. For values of Kg and Kg outside the crosshatched area, no 
sustained oscillations should exist. 

The motion of the roots with variable Kg and Kg are plotted in 
Fig. 1-27, using the data in Table 1-1. The points marked Kg = 0, 

Kg = 0, represent the single-axis LST system solutions. The solid closed 
regions represent the regions of root movements where reliable convergence 
was obtained. The dotted lines include the regions where reliable 
convergence of the iterative scheme was not obtained with 20 iterations, 
and, which can be considered as a boundary area between existence and non- 
existence of solutions. Note that in Fig. 1-27 the two solutions are 
farthest apart when ^ = Kg = Figure 1-27 also shows how the two 



solutions tend to merge just before the solutions disappear. 


ROOT #1 


ROOT n 


h 

K 3 



FREQUENCY 

AMPLITUDE 

FREQUENCY 

AMPLITUDE 

0.1 

0.1 

4.1 

5.07 x TO " 6 

IHPil 

5.99 x io' 7 

0.01 

0.01 

4.26 

5.52 x 10' 6 

mm 

4.58 x 10" 7 

0.001 

0.001 

4.27 

5.56 x 10' 6 

mm 

4.45 x 10" 7 

0.0 

0.0 

4.27 

5.57 x 10“ 6 


4.45 x 10' 7 

0.01 

0.001 

4.27 

5.57 x 10' 6 


4.46 x 10‘ 7 

0.1 

O.OOfv 

4.28 

5.55 x 10' 6 

1.79 

4.59 x 10' 7 

1.0 

0.001 

4.31 

5.24 x 10’ 6 

2.11 

5.79 x 10" 7 

2.0 

0.001 

4.34 

4.91 x 10' 6 

2.43 

7.17 x 10‘ 7 

3.0 

0.001 

4.36 

4.56 x 10" 6 

2.73 

8.69 x 10" 7 

4.0 

0.001 

4.38 

4.20 x 10' 6 

3.01 

1.04 x 10' 6 

5.0 

0.001 

4.37 

3.81 x 10" 6 

3.29 

1.25 x 10' 6 

6.0 

0.001 

4.34 

3.35 x 10" 6 

3.58 

1.53 x 10" 6 

7.0 

0.001 

4.19 

2.6 x 10‘ 6 

3 . 97 ** 

2.06 x io -6 ** 

8.0 

0.001 

3 . 95 ** 

1.50 x 10 " 6 ** 



1.0 

0.001 

4.31 

5.24 x 10' 6 

2.11 

5.79 x 10' 7 

1.0 

0.01 

4.29 

5.19 x 10' 6 

2.23 

5.93 x 10' 7 

1.0 

0.1 

4.13 

4.74 x 10‘ 6 

2.12 

7.71 x 10' 7 

1.0 

0.2 

3.91 

4.18 x 10’ 6 

2.39 

1.03 x 10" 6 

1.0 

0.3 

3.62 

3.47 x 10' 6 

2.62 

1.43 x 10' 6 

1.0 

0.4 

_ 


3 . 2 ** 

2.57 x 10 ' 6 ** 

5.0 

0.001 

4.37 

3.81 x 10‘ 6 

3.29 

1.25 x 10 6 

5.0 

0.01 

4.34 

3.72 x 10" 6 

3.31 

1.30 x 10 6 

5.0 

0.1 

5 . 46 ** 

7.13 x 10 " 6 ** 

3 . 75 ** 

2.25 x 10 6 ** 

1.0 

0.2 

3.91 

4.17 x 10‘ 6 

2.39 

1.04 x 10' 6 

2.0 

0.2 

3.85 

3.64 x 10" 6 

2.79 

1.37 x 10' 6 

3.0 

0.2 

3.69 

2.86 x 10“ 6 

3.28 

1.95 x 10' 6 

0.001 

0.001 

4.27 

5.57 x 10' 6 

1.75 

4.45 x 10' 7 

0.001 

0.01 

4.26 

5.54 x 10' 6 

1.76 

4.56 x 10' 7 
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K 3 

ROOT #1 

ROOT #2 

FREQUENCY 

AMPLITUDE 

FREQUENCY 

AMPLITUDE 

— 

0.1 

4.11 

5.12 x 10' 6 

1.84 

5.79 x IO” 7 

■ 

0.2 

3.92 

4.63 x io' 6 

1.95 

7.52 x IO" 7 

0.001 1 

0.3 

3.7 

4.08 x io' 6 

2.1 

9.82 x IO' 7 

0.001 

0.4 

3.41 

3.42 x IO' 6 

- 

- 

0.001 

0. 5 *' 

2.94** 

2.39 x l(f 6 


- 

0.1 

0.1 

4.1 

5.08 x icf 6 

1.88 

5.98 x io' 7 

0.1 

0.2 t 

3.92 

4.59 x 10" 6 

2.0 

7.79 x IO' 7 

0.1 

0.3 

3.7 

4.03 x io" 6 

2.15 

1.02 x 10' 6 

0.1 

0.4 

3.40 

3.34 x io' 6 

2.38 

1.39 x To’ 6 

0.1 

0.5 

- 

1 

2.82** 

2.22 x IO’ 6 ** 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N^N^G£=0 . 
GfiMMfl= 1.38 00 CD 07 TO= 1.00000D-01 JV- l.OOOOOD 05 

KR= 2.16000D 08 K I = 9 .700000 03 H= 6. OOOOOD 02 

Kl= 1.371 OSD 03 K2= 1 .00000D-01 K3= 1.00000D-01 


JG= 2 . 1 0 0 0 0 D 00 
K0= 5.75835D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .50 0 00000000 

01 

0 .50000000000-05 

1 

0.4733762972D 

01 

0.5433 07 1349D- 05 

2 

0.4526562584D 

01 

0.5454333473D-05 

3 

0.43753230 04 D 

01 

0 .53600155360-05 

4 

0.42709759620 

01 

0 .52599075330-05 

5 

0 . 42032325730 

01 

0. 5 1 83 1873 32 D- 05 

6 

0 .41613238550 

01 

0.51320008560-05 

7 

0.4 13628645 ID 

01 

0.510019 083 1 D- 05 

3 

0.41213130 05D 

01 

0 . 5 08 0665473D- 05 

9 

0.41115734500 

01 

0. 5 0661 8772 OD- 05 

10 

0.41 8522 72880 

01 

0.51180306000-05 

11 

0.41 154884370 

01 

0 . 5 099437398B- 05 

12 

0.41 04483976D 

01 

0. 5 0751 5677 OD- 05 

13 

0.41 08951 0550 

01 

0 .5075 09 0383 D- 05 

14 

0.41 08776583D 

01 

0 .50744499310-05 



ND 32 1 FREQUENCY 3 4.10860D OORAD^SEC AMPLITUDE 3 5.07397D-06 NIT- 14 


LST SYSTEM-NUNER I CRL SOLUTION OF 1+N*61+N*N*6£=0 . 
GAMMA 3 1.380 ODD 07 T0= 1.00000D-01 JV= l.OOOOOD 05 

KP= 2.1600^0 02 K I = 9 .70 GOOD 03 H= 6.Q000QD 02 

Kl- 1.371 OSD 03 K2= 1.00000D-02 K3« 1.00000D-02 


JG= 2 . 1 00 ODD 00 
K0= 5.75835D 03 


ITERATION FREQUENCY 

0 0.41 08597203D 

1 0.4168820145D 

2 0.42 0553 1785D 

3 0 .4226335 066B 

4 0 .42379364720 

5 0.42445 082 OOD 

6 0 .424672 1288D 

7 0 .42532985220 

8 0.42533 l 8862D 

9 0.4253573 1821' 

10 0 .426269652 ID 

11 0. 4255858 0060 

12 0.42552464? ID 


AMPLITUDE 

01 0 .50739665840-05 
01 0 .5248288860D-05 
01 0.5364353532D-05 
01 0.5432259787D— 05 
0 1 0 . 547066358 1 D- 05 
01 0.5492764659D-05 
01 0.5477545295D-05 
01 0.5511 02875 OD- 05 
01 0.551 754 0846B- 05 
01 0.55226 0793 OD- 05 
01 0 .55205888800-05 
01 0 .55124774520-05 
01 0.5516839666D-05 


ND- 1 FREQUENCY 3 4.25522D OORAD^SEC AMPLITUDE 3 5. 521140-06 NIT 3 12 


LST SYSTEM-NUMERICAL S 
GAMMA= 1.38 GOOD 07 TD = 
KF- 2. 16 GOOD 02 KI* 
Kl= 1.37102D 03 K8= 


OLUTIDN OF 1+N*G1+N*N*6£=0. 


i . oooo on- oi 

9.70000D 03 
1 . OOOOOD- 03 


JV= l.OOOOOD 05 
H= 6. OOOOOD 02 
K3= 1.000 OOD- 03 


JG= 2.10000D 00 
K0= 5.75S35D 03 


ITERATION 

0 

1 

2 

3 

4 

5 

6 
7 


FREQUENCY 
0. 42552 17S24D 01 
0.4260676424D 01 
0. 42565348 16D 01 
0.42526 13 07 OD 01 
0.4261846964H 01 
0 .48666934590 0 1 
0 .426953 0544 D 01 
0 .42690428 1 8D 01 


AMPLITUDE 
0.5521 1 42485D-05 
0. 55427958 18D- 05 
0 .53273969681'- 05 
0 . 542 0675622P- 05 
0. 549651 763 OD- 05 
0 .55361 3 5 5 09 D - 0 5 
0 .55578356470- 05 
0 .55615457390-05 


ND= 1 FREQUENCY 3 4.26902B QORAD^SEC HMPLITUBE 3 5.56645D-06 NIT= 7 


Figure 1-12 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G£=0 . 

GAMMA* 1.33000D 07 TO= 1.00 GOOD- 01 .JV= 1 .00 000D 05 JG= £.10 GOOD 00 
KP* £.160000 02 K I = 9.700000 03 H= 6.000000 02 K0= 5.75335D 03 
Kl= 1.37102D 03 K2= 1 . 00000D-Q1 K3= 1.00000H-01 


I TERAT I ON FREQUENCY 


AMPLITUDE 


0 

0.20 0 0 0 0 0 0 0 OD 

01 

1 

0.1 849233 02 4 D 

01 

U- 

0.1 85295 0575 D 

01 

3 

0 . 13671 3923 7D 

01 

4 

0 . 1874 82 954 5 D 

01 

5 

0.1 87889 023 ID 

01 

6 

0 . 188012901 ID 

01 

7 

0 . 1 83 02 0 1 34 6 D 

01 

8 

0.1 88 033435 OD 

01 


0.1 0 0 0 0 0 0 0 0 0 D - 05 
0. 62428486 0211- 06 
0 .5332931 170D- 06 
0.593 06262 14D- 06 
0 .596227 83 06 Ii - 0 6 
0 .59791385860-06 
0 • 5 9 7 9 3 2 8 0 8 4 D - 06 
0.598615391 0D- 06 
0 .59367 05996D- 06 



ND= 1 FREQUENCY- 1 -88042D OORAD/SEC AMPLITUDE* 5.98812D-07 NIT* 3 


LST SYSTEM-NUMERICAL SOLUTION OF 1+M*G1+N^N*G2=0 . 

GAMMA* 1 . 38 0 0 OD 07 T0= 1.0000 OB- 0 1 JV= 1 . 0 0 0 0 OD 05 -JG* £.1000 OB 0 0 
KP= 2.1600 on 02 l - KI* 9. 70 GOOD 03 H* 6.0 000 OD 02 K0= 5.75S35D 03 
Kl= 1.37102D 03 \ K£= 1 . 00000D-02- K3= 1.00000D-02 


ITERATION 


FREQUENCY 


AMPLITUDE 

0 

0. 

1S30421807D 

01 

0.59881 18404D-06 

1 

0. 

18 16575499D 

01 

0.52 04 7436 OSD -06 

d 

0. 

1787791 1 16D 

01 

0 . 4874 09 0332 D- 06 

3 

0 . 

1 774 07148 OD 

01 

0.4720024797D-06 

4 

0. 

176731 0255D 

01 

0 .4645333492D—06 

5 

0. 

176 3 72692 9 D 

01 

0 .46 05378423D— 06 

6 

0. 

17588742 1 7D 

01 

0. 45257934 09D- 06 

“7 

J 

0. 

1758312030D 

01 

0.461 38776 OSD— 06 

3 

0. 

1 758882 166D 

01 

0.45 8 3 0 9 5 3 3 2 D - 0 6 

9 

0. 

1 r' 5 9 6 U 6 o 8 7 D 

01 

0.457761 2 065B- 06 

10 

0. 

175961 7237D 

01 

0.4581 137707D- 06 

1 FREQUENCY* 1 . 759631! 

OORAD^SEC AMPLITUDE* 4.5S129D-07 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +N*N^G2= 0 . 

GAMMA* 1 .33 GOOD 07 TO* 1.00 GOOD- 01 JV= 1.00 GOOD 05 JG= 2.10000D 0 0 
KP* 2.16 0 0 Oil 02 KI= 9. 70 GOOD 03 H = 6. 00 GOOD 02 K0= 5.75S35D 03 
K 1 = 1.371 0£H 03 K2= 1.0 0 0 0 OD- 03 K3= 1.00 0 0 OD- 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 1 759634 0 1 2D 

01 

0 . 458 1 286 7 62D- 06 

1 

0 . 1 754793 059D 

01 

0.451 58 38346 D - 06 

2 

0 . 1 758 0868840 

01 

0.4481 194 31 4D- 06 

o 

0. 174371 6 093D 

01 

0 . 4 3 8 3 3 6 6 5 1 8 D - 0 6 

4 

0.1 7447852 19D 

01 

0.44031761 04D-06 

5 

0.1746624895D 

01 

0 . 4 4 8 3 c. i 0 o otuD U ts 

6 

0.17471 72576D 

01 

0 .4436675832D-06 

7 

0.1 74 7595835 D 

01 

0 . 4 4 4 8 5 7 8 9 8 3 D — 06 

S 

0.17 47828656D 

01 

0 .444676586 1 D— 06 


ND= 1 FREQUENCY* 1 .74799D OORAD^SEC AMPLITUDE* 4.44921P-07 NIT* 8 


Figure 1-13. 
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LST SYSTEM-NUMERICAL SOLUTION OF H-N*61+N*N*G2*0 . 

GAMMA* 1.38 00 OB 07 TO= 1.0 0 000D-01 JV* 1.0 00 GOB 05 JG= 2. 1000 OB 0 0 
KP= 2.1600 OB 02 K I = 9.70 0 0 0 D 03 H* 6.000 0 UD 02 K 0= 5.7 5S35D 03 
Kl= 1.371 Q2B 03 K2= 1.0 GOO OB- 02 K3= 1.0000GB-03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .4269 023 76 SB 

01 

0 .556644825 OB- 05 

1 

0 . 426 043 09 02’B 

01 

0.556 3841 17SB-05 

k. 

0 . 4 2' 6 7 4 3 4 8 6 2 B 

01 

0 .556485936 OB- 05 

3 

0.427236375 IB 

01 

0 . 5564759926D— 05 

4 

0. 42725 12737D 

01 

0.55717 03289B- 05 

5 

0.42727 0092 OB 

01 

0 .557 143743 ID- 05 


ND= 1 FREQUENCY- 4.27286B OORAD/'SEC AMPLITUDE* 5.57162D-06 NIT* 5 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G 1+N^N*G2=0 . 

GAMMA- 1.33 00 OB 07 TD= 1.00 00 OB- 01 JV= 1.00 00 OB 05 JG* 2. 10 00 OB 00 
KP* 2. 1600 OB 02 KI= 9.7000GD 03 H= 6. 00 00 OB 02 K0= 5.75S35B 03 
K 1 = 1 .371 02B 03 ' K2* 1 . 0 0 0 0 OB- 0 1 K3= 1 . 0 0 0 0 OB- 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.427236 04 04B 

01 

0.5571 6 17068D-05 

1 

0 . 4277637363D 

01 

0 . 556 1 946S69B— 05 

2 

0.4275891 1 SOB 

01 

0 .555 0585436 D- 05 

3 

0.4275463586B 

01 

0 .554527 1994D- 05 


NB= 1 FREQUENCY* 4.27574B OORAD/SEC AMPLITUDE* 5.54293D-06 NIT* 3 


LST SYSTEM-NUMERICAL SOLUTION DF 1+N*G1+N*N^G2=0 . 
GAMMA* 1.330 0 OB 07 TO= 1 . 00 0 0 OB- 0 1 JV* 1 .0000 OB 05 

KP* 2. 16 00 OB 02 KI = 9. 7 0 00 OB 03 H* 6. 00 CO OB 02 

K 1= 1.37102D 03 K 2= 1.0 0 000B 0 0 K3= 1 .0 00 0 OB- 03 


JG* 2. 10 00 OB 00 
K0= 5 .7 5835 B 03 


ITERATION 

FREQUENCY 

0 

0 .48757 3 6 6 1 U B 

1 

0 .433 5 23 529 7D 

o 

0.4324421 026 B 

3 

0.431 6383 079B 

4 

0.43 1237298 IB 

5 

0.431 0288669B 

6 

0.43091 1471 OB 

7 

i 

ij . 43 0762225 OB 

R 

0.43 024866 15B 

9 

0 .4307 06 u45 3 B 

10 

0 .43 0693 2 84 SB 


AMPLITUDE 

01 0 .5542934374D— 05 
01 0 . 543 1 663 1 1 2D- 05 
0 1 0 . 5365 104 03 5 B - 05 
01 0 . 53 03 1 6 OS 1 1 D- 05 
01 0 . 5272793 022D- 05 
01 0 . 525743243 OB- 05 
01 0 . 52493 02S55B- 05 
01 0 .5243785475D-05 
01 0 . 52369 04 24 5 D- 05 
01 0 . 52370 03 046D- 05 
01 0. 52363 15323D- 05 


NB= 1 FREQUENCY* 4.30676D OORAB/SEC AMPLITUDE* 5.236S3B-06 NIT* 10 


Figure 1-1 4a . 



32 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N^N*G2=0 . 

GAMMA* 1 . 38 0 0 OH 07 TO* 1.00 0 0 OB- 0 1 -JV= 1.00 0 0 OB 05 Jb= 8.1000 OB 0 0 

KP= 2.16 GOOD 02 KI = 9. 70 00 OB 03 H= 6.QOOOOB 02 K0= 5.75835D 03 

Kl* 1.37102D 03 K2= l.OOOOOB 00 K3= 1 .00 0 0 OB- 03 

HERAT I DN FREQUENCY 

0 0.430 00 OOOOOB 01 

1 0.43 0523 1 336 D 0 1 

2 0.43 0 528 33 5 0 B 0 1 

3 0.4 3 0874487 OB 01 

4 0 . 4 3 0 3 6 3 4 3 5 3 B 01 

NB= 1 FREQUENCY* 4.30865B OORAB/SEC AMPLITUDE* 5.24135D-06 HIT= 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1+nVg1+N*N*G2=0 . 
GAMMA* 1 . 33 0 0 OB 07 TC* 1 . 0 0 0 0 OB- 0 1 JV= 1.00 0 0 OB 05 

K P= 2.16 0 0 OB 02 K 1= 3. 7 000 OB 03 H= 6. 0000 OB 02 

K 1 = 1.37102B 03 K‘2= 2. 0000 OB 00 K3= 1.000 0 OB- 03 

I TERAT I ON FREQUENCY AMPL I TUBE 


0 

0.4 

3 036451 93 D 

01 

0 . 524 1 35 0229B- 05 

1 

0.4 

37596 00 05 D 

01 

0.5 135363 19 OB -05 

2 

0.4 

36 0 0335 0 OB 

01 

0 .5047331 4 12D-05 

3 

0.4 

8 49984778 B 

01 

0 .49764b'653cD— 05 

4 

0.4 

":4 49 053 6 OB 

01 

0.4941 91 3249B-05 

5 

0.4 

342313352B 

01 

Q.4924475958B— 05 

6 

0.4 

3403631 07B 

01 

0 .49 1 52476 05B- 05 

i 

0 .4 

Jo98o4:ju6B 

01 

0 .49 0994755 OB- 05 

3 

0 .4 

3390629 8 3 B 

01 

0 . 4 9 0 6 5 3 6 2 2 3 B - 0 5 


HB= 1 FREQUENCY* 4. 3331 OB QORAB^SEC AMPLITUBE= 4.90415B-G6 NIT= 8 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G2=0 . 

GAMMA* 1 . 33 0 n OB 07 TO = 1.0000 OB- 0 1 J V= 1 . 0 0 0 0 OB 05 

KP= 2.160 0 OB 02 K I = 9.70 0 0 OB 03 H= 6.00 0 0 OB 02 

K 1 = 1.37102B 03 • K2= 3. 00 00 OB 00 K3= 1.0000 OB- 03 

AMPL I TUBE 
0.4904 1521 03B- 05 
0 .4361 0 091 23B- 05 
0.4711 098 35 2D— 05 
0 .46363 0 181 5B- 05 
0 .4599948927B-05 
0.4 53 1 578 3 9 4 B- 05 
0.4571 77352 1 B- 05 
0.45661331 05B- 05 
0. 45625484 17D- 05 

NB* 1 FREQUENCY* 4.36199B OORAD^SEC AMPLITUDE* 4.55997D-06 NIT= 8 


ITERATION FREQUENCY 


0 

0 .43831 00273D 

01 

1 

U .44 079 4 7 3 3 0 B 

01 

2 

0 .433733349 OB 

01 


0.4375 5 0 6 9 4 9 D 

01 

4 

0. 4369644 072B 

01 

cr 

■j 

0 .43666 1333 IB 

01 

6* 

0 .43649 1822 5 D 

01 

*7 

i 

0.4 36381 55 07B 

01 

o 

U . 4 36 893 9 4 y ID 

01 


•JG= 2.1 00 0 OB 0 0 
K0= 5.75S35B 03 


JG= 2. 10 00 OB 00 
KG* 5.75S35B 03 


AMPLITUDE 
0.520000000 OB- 05 
0.52233 0390 IB- 05 
0 . 582245 3439B- 05 
0.524574872SB— 05 
0.5243241 015D- 05 


‘tSSSS^ 


Figure 1-145. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N^N^G2=0 . 
bftMMR= 1 .38000H 07 TO= 1. 000000-01 JV= 1.00000D 05 

KP= £.16 00 0 D 02 K I = 9.70 0 0 OD 03 H= 6.000 0 00 02 

Kl= 1 .3710211 03 K2= 4.00000D 00 K3= 1 . OOOOuB-03 


JG= 2.1000 OE 00 
K0= 5.758350 03 


ITERATION 

0 

* 1 


FREQUENCY 
0.436 1988 1750 01 
0. 443244 0367D 01 
0.44056315540 01 
0.439 1262 01 60 01 
0.43842438520 01 
0 . 438 06287320 01 
0.43 7859 033 OD 01 
0.43773946740 01 
0.43766507050 01 


AMPLITUDE 
0 .45599656840-05 
0 . 452593 06 160- 05 
0 .436 0942 0230- 05 
0.428 0822 0260-05 
0 .42421 089640-05 
0.42224945990-05 
0.421 19239210-05 
0 . 42 053939390- 05 
0.42 022 1326 ID- 05 






ND= 1 FREQUENCY = 4.37615D OORAD'SEC AMPLITUDE** 4.199330-06 


N I T= 


LST SYS T E M-NUME P I CRL SOLUTION OF 1+N*G 1 +N^N*G2= 0 . 


GHMMR= 1 .38 

0 0 OD 

07 TO** 1.000000- 

01 JV= 1.000000 05 

JG= 

2.100000 

KP= £.160000 02 

K I = 9.700000 

03 H= 6.000000 02 

K0= 

5.758350 

Kl= 1.37102 

0 03 

\ K2= 5.000000 

00 K3= 1.000000-03 



ITERRT ION 

FREQUENCY 

AMPLITUDE 



0 

0 . 4-876 1486130 01 

0 .41 9933 0 0270- 05 



1 

0.44462636900 01 

0.41725291280-05 



2 

0.44107233930 01 

0 . 39862393760- 05 



■Hk 

0.4 

39231 6296 D 01 

0.38973123510-05 



4 

0.4 

3334979140 01 

0.38548266220-05 



5 

0.4 

3789680840 Ul 

0.3833214 1800-05 



6 

0.4 

3763990410 01 

0 . 382 1 3926590- 05 



7 

0.4 

3748983080 01 

0 . 3 3 1 4 5 3 6 9 040- 0 5 



o 

0.4 

3740363500 01 

0.38103502530-05 



9 

0.4 

3735670320 01 

0.38076891610-05 



NO= 1 FREQUENCY- 4.373330 OORRO 

/SEC RMPLITIJHE= 3.805 

950-06 

N I T= 9 

LST SYSTEM- 

NUMERICAL SOLUTION OF 

1+N*G1+N*-N*G£=0 . 



GAMMA** 1.38 

00 00 

07 TO = 1.000000- 

01 JV= 1 . uOOOOD 05 

•JG= 

2 . 1 00000 

KP= 2.160000 02 

KI= 9.700000 

03 H= 6.000000 02 

K0= 

5 . 758350 

K 1= 1.371020 03 

K2= 6. 0 0 GOOD 

00 K3- 1.000000-03 



ITERATION 

FREQUENCY 

AMPLITUDE 



0 

0.4 

3733253690 01 

0 . 38 0595 05940- 05 



1 

0 .444 3r. x y ij y : *: 0 01 

0 .37862256920- 05 



o 

0.4 

3937160140 01 

0.356394471 10-05 



3 

0 . 4 

-;kk 9 085630 Ul 

0.34560153840-05 



4 

0 .4 

3545475270 01 

0 . 34 05764 0 700- 05 



5 

0 .4 

34S162746D 01 

0 . 3 3 7 9 9 1 6 3 6 5 D - 0 5 



6 

0.4 

3444856440 01 

0.33653822720-05 



7 

0.4 

3422667430 01 

0.3356649 1440- 05 



y 

0.4 

3411479330 01 

0 .33513650170-05 



9 

0.4 

3408331 140 01 

0 • 3348264 1 c'SD— 05 




ij i.i 
03 


ND= 1 FREQUENCY= 4.340810 OORAD/SEC AMPLITUDE** 3.346470-06 NIT= 9 


Figure -14c 
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L3T SYSTEM-NUMER I CAL S 
GAMMA* 1 . 33 0 0 00 07 TD= 
KP= 3.1600 OD 03 KI = 
Kl= 1.371 OSD 03 Kc- 


□LUTIDM QF 1+N^G1+N*N*G2*0. 

1 . 0 0 0 0 0 D - 0 1 J V = 1.0 0 0 0 0 D 0 5 

9.70 0.0 OD 03 H= 6. 0000 OH 03 

7.0 0 0 0 OD 00 K3= 1 .00 GOOD- 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .434 03 1233 ID 

01 

0 . 3-34 64 7 -3 -3 73 0— U 5 

1 

0 .44 08645 04 9D 

01 

0 . 3 o 2 5 7 8 975 8 D — 0 5 

£ 

0 .43250391 46D 

01 

0.30135972570-05 

. p m 

0. 42627 01797D 

01 

0. £323294 047D- 05 

4 

0 .4233 1322 OSD 

01 

0 .2738939 036D- 05 

5 

0.4216530616D 

01 

0 . 2683925 0650- 05 

6 

0.42 049353 06D 

01 

0. 2655545 194D- 05 

( 

0.41 9 7396293 D 

01 

0.26341 0 01 5 OD- 05 

o 

0.41 3375 1327D 

01 

0 .26 1 26 1 9774D- 05 

9 

0 .4 139724 044 D 

01 

0.26 0644 1657D-05 

10 

0.419121 1536H 

01 

0 .26 0444645 OD— 05 


JG= 2.100 0 OD 00 
K0= 5 . 75 3 35 D 03 


ND= 1 FREQUENCY* 4.191370 OORAD-'SEC AMPLITUDE* S.60430D-06 MIT* 10 


LST S 
gamma 
KP* 3 
Kl* 1 


YSTEM-NUMER I CAL 3 
= 1 .33 OOOrt i 07 TD= 
.16 0 0 OB 02 KI = 

.371 OSD 03 K2= 


□LUTIOM UF 1+M*G1+N^M*GS=0. 

l.oo o u OD- o i J v= i .no o o on 05 

9.700 0 0 D 0 3 H = 6.0 0 0 0 0 D 0 2 

S.OOOOOD 00 K3* 1.0G00GD-03 


JG= £. 1 000 OD 00 
K0= 5.75835D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.4191 366 58 6 D 

01 

0 . 26 046: n 0965B— 05 

1 

0 .432083 0567D 

01 

0.275544992011-05 

L- 

0 . 3 8 85654674 D 

01 

0.1 5372275460- 05 

3 

0 « 3 9 9 8 6 8 6 9641) 

01 

0.1 777 0982 03D- 05 

4 

0.410121 1254D 

01 

0.20181 97669B- 05 

5 

0 .4 049 0681 5 7D 

01 

0. 21 337527 15D- 05 

6 

0.421 865 0868D 

01 

0 .84488 08 0 1 60—05 

1 

0.5328541 039D 

01 

0.51215428030-05 

8 

0 .4932374 097D 

01 

n _ 4 o '-i 1 1 9 '-t 6, 1 6. i n — ■ 1 1 s 

9 

0 .4580742990H 

01 

0 .8898 763479 D— 05 

10 

0.42726331 17B 

01 

0.2 8 3 8 9 3 0 8 5 9 D - 0 5 

11 

0 .4 726 01 27 GOD 

01 

0.3527236 01 3D- 05 

12 

U .464846203 8 D 

01 

0 .22062 1 42S7B- 05 

13 

0.41 02385 16 ID 

01 

0 .171 26 07528D- 05 

14 

0.4511 00094 OD 

01 

0 . 8 8 0 1 ti 5 1*' 1 0 1 D — U5 

15 

0 .4 185 049431 D 

01 

0 • 8 1 4 39 9 5 8 4 0 D U 5 

16 

0.433792 154 ID 

01 

0 . 2566385 783D- 05 

17 

0.45847 1929 ID 

01 

0 . 3 1 72396 064D- 05 

18 

0.441 031 31 36D 

01 

0.25954 176 160-05 

19 

0 .39 01 53372 OD 

01 

0.14570902020-05 

2 0 

0 .3951355 3 3 9D 

01 

0.15 0 0 0364790- 05 


MD* 1 FREQUEMCY* 3.951360 OGRRD-'SEC AMPLITUDE* 1.500040-06 MIT* £0 


Figure 1-1 4d 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N^G 1+N<*N*G2® 0 . 

GAMMA® 1.3800011 07 TO® 1 . 0 0 00 OB- 01 JV= 1.00 GOOD 05 . JG= 8.1 0 00 OD 00 

KP= £.1600011 02 KI = 9.70 0 000 03 H® 6.0000GD 02 K0= 5.753350 03 

Kl = 1.371020 03 K2= 1.000000-02 K3= 1.000000-03 


ITERATION 

0 

FREQUENCY 
0 . 174 7994 349D 

01 

AMPLITUDE 
0.444921 0692D-06 

1 

0.1 75 097 0353D 

01 

0 . 445962 8 3 5 6 D - 06 

* 2 

0.1 75 1 578924D 

01 

0 . 44586 1 u5c!4D— Ub 

3 

0 . 1751663355D 

01 

0 .4463 09 084 OD- 06 

4 

0 . 1751762687D 

01 

0 .4463992 173D- 06 


ND = 1 FREQUENCY® 1 .751831) OORAB/SEC AMPLITUDE® 4.46504D-07 NIT® 4 


• LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N^N^6£=0 . 

GAMMA® 1.33 00 OD 07 TO® 1.00000D-01 JV= 1.0 000 OH 05 

KP= 2. 160 ODD 02 KI= 9.700 ODD 03 H= 6.0 000 OD 02 

Kl= 1.37102D 03 K£= 1.00000D-01 K3= 1.0 000 OB- 03 

AMPLITUDE 
0.4465041 105D-06 
0.451 01 01543D-06 
0 .4546263 062B- 06 
0 .456764443 1 D— 06 
0 . 457744 7 044 D- 06 
0.4579707749D-06 

ND= 1 FREQUENCY® 1.73701D OGRAD'SEC AMPLITUDE® 4.53193D-07 NIT® 5 


ITERATION 

FREQUENCY 


0 

0 . 1 75 v l 334 1 76D 

01 

1 

0 . 176^294 195D 

01 

o 

I— i 

0 . 1 7789 04699D 

01 

3 

0 . 1 734 00 9536 D 

01 

4 

0.1 7866667 1 2D 

01 

5 

0 . 1 7S6379223D 

01 


J6= 2.10000B 00 
K 0= 5 • 75835B 03 


LST SYSTEM-NUMERICAL S 
GAMMA® 1.33 GOOD 07 TO® 
KP= 8.1600 OD 02 KI = 
K 1 = 1.37102D 03 K2® 


□LUTION OF 1 +N+61 +N*N+G2= 0 . 
1.0000 OD- 0 1 JV= 1.00 0 0 OD 05 

9.70 GOOD 03 H= 6.00000D 02 

1 . 0 0 0 0 OD 0 0 K3= 1.0000 OD- 03 


J6= £.10 GOOD 00 
K 0= 5 . 75S35D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.173701 001 2D 

01 

0 .458 1926 03 ID- 06 

1 

U . 13664 6 3 1 8 8 D 

01 

0.331451 1 173D- 06 

c! 

0 . 19539991 2 OD 

01 

0.41 4834897 OD— 06 

3 

0 . 8 02 1 29 0853D 

01 

0 . 4S25286553D- 06 

4 

0.2 0645237 09D 

01 

0 . 5264393 0 04D- 06 

5 

0.20891 01 053D 

01 

0.551 4356 079D- 06 

6 

0 .21 01973469D 

01 

0 . 56443 0 03 05B- 06 

7 

0 . 8 1 0326652 1 D 

01 

0.57 0669 02 02D- 06 

£j 

0.211 0689 09 AD 

01 

0 .57351 03622D-06 

9 

0.2111 477387D 

01 

0 . 5755678257D- 06 

10 

0.21121 35 359 D 

01 

0 .57648543 7 2D- 06 

11 

0 .21 1249181 3D 

01 

0 . 57721 922 62 B- 06 

12 

0.2 11 27688 19D 

01 

0 . 5 7 7 64 324 9 4 D— 06 


ND= 1 FREQUENCY® 2.11296D OORAD'SEC AMPLITUDE® 5.77976D-07 NIT® 12 


Figure l-15a. 



LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*Gl+ N*N*G2=G . 


GAMMA* I • 3800 CD C7 tq = UOOOOOD-OL J V = l.QOQUQD 05 JG= 2,100000 00 


v? MCI IM - i • jou 

KP= 2.160000 

02 K I = 9. 

700CGD 

C3 H* 6 • QOOOOD 02 

K0= 5*7 5 835 D 03 


03 K 2= L. 

000000 

00 K 3= 1.000 00D-03 


iteration 

FREQUENCY 


AMPLITUDE 


0 

0.2100000000D 

Oi 

0.5 8000000000-06 


1 

C. 21058192980 

01 

0.5 76440 82180-06 


2 

0.21 106834450 

01 

0.57779394370-06 


3 

0.21 13400637D 

01 

0. 57848853980-06 


• 4 

0.21 134585760 

01 

0.57875519580- 06 



NO* I FREQUENCY* 2.11352D OORAD/SEC AMPL I TUDE = 5. 78829D-07 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF i+N*Gl+N*N* G2 = 0 . 

GAMMA= 1,380000 07 TO* l.OOQOOP-Ol J V= l.QOQQOP 05 JG* 2.1000QD QO 

KP = 2. 160000 C2 K I * 9.7COOOD 03 H= 6.000000 02 K0= 5.758350 03 

Kl = 1.371020 03 K 2= 2.00000D 00 K 3= 1.00000D-03 


ITERATION 

0 

i 

FREQUENCY 
0. 21 1352 23 68D 
0.21987837660 

01 

01 

AMPLI TUDE 
0.5788 29 4272D- 06 
0. 47457342560-06 


2 

0.22862730950 

01 

G. 55415053110-06 


3 

0.23516917170 

01 

0.62337062680-06 


4 

0. 2 3 9234 06 82D 

01 

0.66771529620-06 


5 

0.24145437180 

01 

0.6 9194946710-06 


6 

0.2425574373D 

01 

C. 7037403 84 8D- 06 


7 

0.243C351540C 

01 

0.7C88785615 0-06 


8 

0.24 31 88 08 54D 

01 

C. 7 1202441000-06 


9 

0.24328 39304D 

01 

C . 7 1389960070-06 


10 

0. 24334826530 

01 

C.715106C090D-06 


11 

0 # 2 4 3 39 3 54 86D 

01 

C.7 1593866780-06 


12 

0.24 342 70 5 3 8D 

01 

0.71653880260-06 



ND = 1 FREQUENCY* 2.434530 OORAD/SEC AMPLITUDE* 7. 16991D-07 NIT* 12 


LST SYSTEM-NUMERICAL SO LUT I C N CF 
GAMMA*. 1.380000 07 TO* i.OOOOOD- 
KP = 2.160000 02 K I = 9.70000D 

l + f\*Gl + N*N*G2*0. 

01 J V* l.OOOOOD 05 

03 H= 6.00000D 02 

JG* 2.10000D 
KO* 5.758350 

00 

03 

K 1= 1,371020 

03 K 2 = 3 . < 

OOOOOC 

00 K 3 = 1.00000D-03 



ITERATION 

FREQUENCY 


AMPLITUDE 



0 

0.24345273020 

01 

C. 7 1699083600-06 



1 

0.25307037250 

01 

0.675391 79420-06 



2 

0.26 I5049668D 

01 

C. 74730169010-06 



3 

0.267181895 70 

01 

C. 80575802490-06 



4 

0.2703849298D 

01 

0.839567 19440-06 



5 

0. 27 198 19820D 

01 

C .85575344710-06 



6 

0.2 7266 5605 8D 

01 

C. 86176561730-06 



7 

0.2728L7 34610 

01 

0.8639 00 6 04 90-06 



8 

0.272891911 70 

01 

C. 86589763960-06 



9 

0.272960325 ID 

01 

0. 8669499962D-06 



10 

0.27 3004 545 ID 

01 

0.8 679 340734 D- 06 



11 

0.27304218240 

01 

C .8 6858662080-06 




ND= 1 FREQUENCY* 2.73071D QOPAO/SEC AMPLITUDE* 8.691500-07 NIT= ll 


Figure 1-1 5b. 


37 


1ST SYSTEM-NUMERICAL solution of 
GAMMA = I • 300 COO C7 TO* l.OOOQOD- 
KP= 2.I6000D 02 KI= 9.70GQ0D 

1 + N*G1+N*N*G2=0 . 

01 JV= 1 . 000 ODD 05 

03 H = 6.00000D 02 

JG= 2.100000 
K 0= 5.75835D 

00 

03 


KI= L • 371023 

03 K 2= 4.' 

QQOGOD 

00 K 3= 1.00000D-03 




ITERATION 

FREQUENCY 


amplitude 




0 

0*27 307 1 12410 

01 

0.8691 49 8890 D- 06 




1 

0. 283 276 1272D 

01 

C .869056 2992D-06 




2 

0.29 140723690 

01 

0.940441 74470-06 





3 

0. 29652075090 

01 

0.993123 0333D-06 




4 

0. 29 927995 95Q 

01 

0*1021781 45 2D- 05 




5 

0 • 30 06 306 5 16D 

01 

0. 1C34869396D-05 





6 

0*3011 1433820 

01 

0.1G3300193BD-05 




7 

0 • 30 1 160 5 30 6D 

01 

0. 1 C40285461Q-05 




8 

0.301229^6140 

01 

0. 1040916602D-05 





9 

0*30 126669 6UD 

01 

0.1041986665 D-05 




10 

0*30 13.04 1 149D 

01 

0.1 C4258 96 02D-05 





ND= 1 FREQUENCY* 3.01333D OORAD/SEC AMPLITUDE* 1.04321D-06 NIT= 10 


LST SYSTEM-NUMERICAL SOLUTION OF 
GAMMA* 1.38CC0D C7 TO* l.OOOOOD- 
KP= 2. 16000D 02 KI= 9.70000D 

1 + N*G1 + N*N*G2 = 0 . 

01 J V* l.OOOOOD 05 

03 H= 6. OOOOOD 02 

JG= 2.10000D 
K 0= 5 .75 83 50 

00 

•03 

K 1= 1. 37102D 

03 K 2= 5.1 

300 OOD 

00 K3* 1 .000000-03 



ITERATION 

FREQUENCY 


AMPLITUDE 



0 

0.30 13329729D 

01 

0.1 043207423D-05 



1 

0.31 158 8 594 8D 

Cl 

0.1 C65 99 22 250-05 



2 

C. 31957166540 

01 

0.1 1425688390-05 



3 

0 . 3245957 530C 

01 

0.1 1987193770-05 



4 

0 .32 7357447 ID 

01 

0. 123019 1890D-05 



5 

0.32 86982993D 

01 

0.1 2441944220-05 



6 

0. 32899643840 

01 

0.12456812450-05 



7 

0.32 90445 5 500 

01 

0.1 24873 1532D-05 



8 

0 • 3 29 1 33 3009D 

01 

0. 1 2493769450-05 



9 

0.3 29 17659 32 D 

01 

0.12508261720-05 



10 

0. 329224476 7D 

01 

0. 12515835020-05 




ND= 1 FREQUENCY = 3.29261D OORAD/SEC AMPLITUDE* 1.252440-06 M1T= 10 


Figure 1 -15c. 
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L ST SYSTEM- NUMER ICAL SOLUTION OF 
GAMMA= 1.38CCCD 07 TQ = l.OOOOOD- 
KP= 2.160000 02 K I - 9.70000D 

l+N*Gl+N*N*G2=0 . 

01 J V- 1 *000000 05 

03 H = 6.00000D 02 

JG= 2.10C00D 
K0= 5.750350 

00 

03 

K L- 1.371020 03 K2= 6. 

COOCOO 

00 K3= 1.000000-03 



ITERATION 

FREQUENCY 


AMPLI TUOE 



0 

0.32926089030 

01 

0.12524417490-05 



1 

0.33931182040 

01 

0.12880153440-05 



2 

0. 34766480C4D 

01 

0.1382 1226540-05 



' 3 

0.35324191430 

0 1 

0*14555935280-05 



4 

0.35645857530 

01 

0.14992048020-05 



5 

0.35782924460 

01 

0. 15150643390-05 



6 

0.357995951 7Q 

01 

0.151789 12160-05 



7 

0.3581068 133D 

01 

0. 1520970971 D- 05 



8 

0.358215 34120 

01 

0. 1 5230940250-05 



9 

0.35829991 13D 

01 

0.15251788380-05 



10 

0.35 837765440 

01 

0.15268278690-05 



11 

0.35844360030 

01 

0. 1528349560D-05 




ND= 1 FREQ U 

ENC Y= 3.585G3D CORAD/SEC AMPLITUDE = 1. 529650-06 NI T= 11 


L S T SYSTEM-NUMERICAL SOLUTION OF 

1+N*G1 + N*N*G2 = 0. 


GAMMA- 1.38 GOOD 07 T 0= 1 .1 

30000D- 

01. JV= 1 .000 ODD 05 JG= 2. 10000D 

00 

K P = 2 . 1 60 COD 

02 K I - 9.' 

70000D 

03 H = 6.000000 02 K0 = 5.758350 

03 

K 1= 1.371020 

03 K2= 7.1 

OOOGOD 

00 K 3= 1.00000D-03 


ITER A t IC-N 

FREQUENCY 


AMPLI JUDE 


0 

0. 3 5 8502945 50 

01 

0.152S65C116D-05 


1 

0.36835499180 

01 

0.1 575826533 D- 05 


2 

0.37823023170 

01 

0.1 7157825550-05 


3 

0.38603849890 

01 

0.18458471720-05 


4 

0.39 1008 43 6 30 

01 

0.193061 7649D-05 


5 

0.39215573620 

01 

0.19473093000-05 


6 

0.39279550100 

0 1 

0. 15590148480-05 


7 

0. 3931729440D 

01 

0. 19900152340-05 


8 

0.39497523190 

01 

0.2C1402 1481D-05 


9 

0.39509208140 

01 

0. 20242424210-05 


10 

0. 39547837670 

01 

0 • 2 02 743 38 04 D- 05 


11 

0.39560248320 

01 

0.2C3302 1197D-05 


■ 12 

0.39578890570 

01 

0. 2035 75 73530-05 


13 

0.39591420550 

01 

0.2039619541 D- 05 


14 

0.39605362450 

01 

0.2C4246 59 0 60- 05 


15 

0.396 17 190700 

01 

0.2C45558664D-05 


16 

0.39628987570 

Cl 

G.2C48246234D-05 


17 

0.39639796620 

01 

0.2C509056 35Q-0 5 


18 

0.39650208010 

01 

0.20533612810-05 


19 

0.39660004090 

01 

0.2C55723118Q-05 


20 

0.39669361610 

01 

0.2057952294D-05 



ND= 1 FRFOUENCY* 3.96694D OORAQ/SEC AMPLITUD6= 2.05795D-06 NIT= 20 


Figure 1-15d. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N^G1 +N*N^G3=0 . 


GAMMA® 1.38 GOOD 

07 TO® 

1 . 00 GOOD- 01 

JV® 

1.00 GOOD 05 

JG= £.10 GOOD 

KP= 3. 160 ODD 03 

K I = 

9.7000 OD 03 

H= 

6.00 0 0 OD 03 

K 0® 5.75835D 

Kl= 1.37103B 03 

K3= 

1.00000D 00 

K3= 

3. oooo on- oi 



ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 39 0 0 0 0 0 0 0 OD 

01 

0.430000000 OD- 05 

1 

0.391 1764448D 

01 

0.41 98 1 5436 9D- 05 

3 

0.391 0445465D 

01 

0.41 856 1 0396D- 05 

■3 

0 . 39 07388 863 D 

01 

0.41 75505 0 1 7D- 05 

4 

0.39 097 15946B 

01 

0.4174333359D-05 


ND= 1 FREQUENCY® 3.91 0 04B OORAD^SEC AMPLITUDE® 4.17373D-06 NIT® 4 


LST SYSTEM-NUMERICAL SOLUTAON OF 1 +N*G1 +N*N^G3=0 . 

GAMMA® 1.33 GOOD 07 TO= 1.00000D-01 JV® 1 .0 000 OD 05 JG= 3.1 0 0 0 OD 00 
KP= 3.1600 OD 03 KI= 3.70 GOOD 03 H= 6. 00 GOOD 03 K0= 5.75835H 03 
Kl = 1 .37103H 03. K£= 3.00000D 00 K3= 2.00000D-01 


ITERATION 


FREQUENCY 


AMPLITUDE 

0 

0. 

391 0 03 8835 D 

01 

0.4 17378 1433D-05 

1 

0 . 

394485861 3D 

01 

0.4 03439 1937D- 05 

£ 

0. 

3*?0£493997D 

01 

0 .33345S9333B-05 

3 

0. 

388 0 0 1 539 7 D 

01 

0.373844 5 8 3 9 D - 0 5 

4 

0. 

386933 556 7 D 

01 

0.369395 1 635D- 05 

e- 

0. 

8 a 6o9 6 3 8 U 1 D 

01 

0 . 367 0 1 44395D— 05 

b 

0. 

3861 01S733D 

01 

0 .-365 •' c* 6 5 6 6 5 D — U 5 

, —9 

1 

0. 

3859879 764D 

01 

0 .365 078 1946 D- 05 

8 

0. 

3858194339 D 

01 

0 .36464 1 0353D-05 

9 

0. 

3857484 154D 

01 

0 .3643544834D- 05 


ND= 1 FREQUENCY® 3.85700D OORAD^SEC AMPLITUDE® 3.64157D-06 NIT= 9 
LST SYSTEM-NUMERICAL SOLUTION OF 1 +N^G1 +N+N+G3® 0 . 


GAMMA® 1.38 GOOD 

07 TD= 

1 . OOOOOD- 

01 

JV= 1.0 00 OOD 05 

JG= 3.10 0 0 OD 00 

KP= 3.160 OOD 03 

K I = 

9.70000D 

03 

H= 6. 00 GOOD 03 

K0= 5 . 75835D 03 

K’l = 1 .37103D 03 

KS= 

3 . oo oo on 

00 

K3= 3.0000 OD- 0 1 



ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 3 3 5 6 9 9 7 8 8 3 D 

01 

ij m 864 1 5 6 9 8 ij 8D — Li 5 

1 

0.33681 51 877D 

01 

0 . 344635645 1 D- 05 

3 

0 .3795 14 01 39D 

01 

0 . 3 1 783 07333D- 05 

. 3 

0 .374474385 OD 

01 

0.3 03338581 3D- 05 

4 

0 .373 0805 047D 

01 

0 .395 0 6 9 4 0 3 3 D - 0 5 

5 

0.37 085365 18D 

01 

0.3914 3 9 9 7 0 3 D - 0 5 

6 

0 . 37 0 0756579D 

01 

0 .339191346 3D -05 

7 «J 

0 . 3 6 9 6 6 9 5 8 7 D 

01 

0 .38793 0881 ID- 05 

3 

0.368941 035 ID 

01 

0 .8864645745D-05 

9 

0 .36934 4 983 9D 

01 

0.38639051 16D- 05 


ND= 1 FREQUENCY® 3.693C4B OORAB/SEC AMPLITUDE® 3.86133D-06 NIT® 9 


Figure 1-16. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N^N^G£=0 . 

GAMMA* 1 .380000 0? TO* 1.00000D-01 -JV= 1.000000 05 JG= 2.10000D 00 
KP= 2.160000 02 KI= 9.700000 03 H= 6.000000 02 K0= 5.75835D 03 
Kl = 1.371020 03 K2- 1.00000D 00 K3= 2.000000-01 


ITERATION 

FREQUENCY 


AMPLITUOE 

0 

0 .24000000000 

01 

0. 1 0 0 0 0 0 0 0 0 OD- 05 

1 

0 .240947944 1 D 

01 

0.10467125640-05 

2 

0.23992907240 

01 

0. 10338041690-05 

3 

0.2394536278D 

01 

0.10357253560-05 

4 

0.2390864967D 

01 

0.10313689430-05 

5 

0 .239507685 00 

01 

0.1 046663 0 1 70—05 

6 

0.23984041350 

01 

0.10513670190-05 

7 

024 0 05393260 

01 

0.1 0424853980-05 

8 

0 .23925323460 

01 

0.10313138350-05 

9 

0.23929046130 

01 

0.10373983360-05 

10 

0.23934011180 

01 

0.10406744110-05 

11 

0.23936038340 

01 

0.1 04072 0365D- 05 


NO* 1 FREQUENCY- 2.39389D 0 ORAL/SEC AMPLITUDE* 1.041350-06 NIT= 11 


LST S YSTEM-UUMEP I CflL SOLUTION OF 1 +N*Gl+HoN*G2=0 . 

GAMMA* 1 .330000 07 T0= 1.000000-01 JV= 1.000000 05 JG= 2.100000 00 

KP* 2.16000ft 02 KI= 9.700000 03 H= 6.00000D 02 K0= 5.758350 03 

Kl= 1.371080 03 K2- 2.000000 00 K3* 2.000000-01 

ITERRTION FREQUENCY AMPLITUDE 


0 

0.23933361850 

01 

0.1 041351 79 ID— 05 

1 

0. "5455332780 

01 

0.11 151902020-05 

2 

0.26589495670 

01 

0.12224496680-05 

3 

0.27300435520 

01 

0 . 12988748 120-05 

4 

0.27687369 050 

01 

0.13413889240-05 

5 

0.27876453620 

01 

0 . 1 36 12 0737 OD- 05 

6 

0.27945522190 

01 

0 . 13666908380—05 

7 

0.27953S2620D 

01 

0.13687160080-05 

8 

0.27961869060 

01 

0.1370021 1 050-05 


ND= 1 FREQUENCY* 2.79679D OQRAD/SEC AMPLITUDE* 1.371270-06 NIT® 3 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G2* 0 . 

GAMMA* 1 .33 00 00 07 TO* 1.00 GOOD- 01 JV= 1 .00 00 0D 05 -JG= 2.100 0 00 00 
KP* 2.160000 02 KI= 9.700000 03 H= 6.000000 02 K0= 5.758350 03 
Kl= 1.371020 03 K2= 3.0Q000D 00 K3= 2.000000-01 


ITERATION 

FREQUENCY 


AMPLITUOE 

0 

0.27967859540 

01 

0.13712719940-05 

1 

0.29573095321) 

01 

0 . 15029147890-05 

£ 

0 . 3 082893 02 00 

01 

0.16542507510-05 

3 

0.31724009810 

01 

0.17798663200-05 

4 

0.322981 09 08D 

01 

0 . 18654885140-05 

5 

0.32552671 09 D 

01 

0.19008 0 37390- 05 

6 

0 .32604438330 

01 

0.19088475590-05 

7 

0 . 3264 0 1 8726D 

01 

0.19167072710-05 

8 

0 .32667987780 

01 

0.1 92 35 097800-05 

9 

0 .326928641 ID 

01 

0.19287213530-05 

10 

0.327 137 12 76D 

01 

0.19333620760-05 

1 1 

0 . 327323 0327D 

01 

0 . 1 9373 063840- 05 

12 

0.32748559800 

01 

0 . 194 082676 OD- 05 

14 

0.32788208080 

01 

0.19489235980-05 

15 

0.32788 0272 ID 

01 

0.19492343150-05 

16 

0.32798762540 

01 

0.19515184270-05 

1? 

0 .32308577690 

01 

0.19535996010-05 

13 

0.32817582330 

01 

0.19555080650-05 


ND= 1 FREQUENCY* 3.23259D OORAD^SEC AMPLITUDE* 1. 957260-06 NIT* 18 


Figure 1-17 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G2*0 . 

GAMMA* 1 .38 GOOD 07 TO= 1 . 00000D-01 JV= 1 . 00000D 05 -Jb* 2. 1 00 0 Oil 0 0 
KP= £.16 0 000 0£ KI= 9.70 0 0 OB 03 H= 6.00000D 02 K0= 5.75835D 03 
Kl* 1.3710£Ii 03 K2= 1.00000D-03 K3= 1.00000D-03 


ITERATION 

FREQUENCY 


0 

0 .43 0000 0 000 D 

01 

1 

0.4271 0S6058P 

01 

2 

0 .4265786694 D 

01 

3 

0. 426875 0642D 

01 

4 

0.429 0753 07 ID 

01 

5 

0.4278 01546311 

01 

6 

0 .4271 0731 96D 

01 

7 

0.4272476170D 

01 


AMPLITUDE 
0 . 55 0 0 0 0000 OD- 05 
0 . 55 08 161 796II- 05 
0 . 5527424953 0- 05 
0.555219758911-05 
0.5645724 07811-05 
0 . 56 04 1 79382D- 05 
0. 55324804 180-05 
0. 5579216 017D-05 



ND= 1 FREQUENCY- 4.27381D OORRD/SEC AMPLITUDE* 5. 5788411-06 NIT* 7 


LST SYSTEM-NUMERICAL S 
GAMMA* 1 .380000 ’i07 TO* 
KP= £ . 1 6 0 0 OD 02 K I = 
K 1 - 1.371 OSD 03 K2= 


□LUTION OF 1+N*G1+N^N^G2=0 . 

1.00 GOOD- 01 s JV= 1 .00 0 0 OD 05 

9.7000 OD 03 H= 6.0 0 0 0 OD 02 

1.0000 OD- 03 k‘3= 1 . 0 0 GOOD- 02 


JG= 2.1000 OD 00 
K0= 5.75835D 03 


I TERAT I ON FREQUENCY 

0 0.427331 1433D 01 

1 0.42651 6528 OD 01 

2 0.425689 19 OSD 01 

3 0.425857 191 5D 01 


AMPLITUDE 
0 .5578838989D-05 
0 . 556 1 562337D- 05 
0 . 553 98 5 4 734 Ii- 05 
0 . 5537666565D-05 


ND= 1 FREQUENCY* 4.25947D OORAD^SEC AMPLITUDE* 5.53740D-06 NIT* 3 


LST SYSTEM-NUMERICAL S 
GAMMA* 1 .3800 01' 07 TO* 
K P = 2.16 0 0 0 D 0 2 K I = 

Kl* 1.371 OSD 03 K2= 


□LUTION OF 1 +N+G 1 +N+N+G2* 0 . 

1 . 0 0 0 0 0 n— 0 1 J V = 1 . 0 0 0 0 0 D O 5 

9 .700 0 OD 03 H= 6 . 0 0 0 0 OD 02 

1 . 0 0 0 0 OD- 03 K3= 1 . 0 0 0 0 OD- 0 1 


Jb= 2.10000D 00 
KG* 5.75835D 03 


ITERATION 

0 

FREQUENCY 
0.42594 7337 ID 

01 

AMPLITUDE 

0. 5537397 133D- 05 

1 

0.41 962 12759D 

01 

0 . 5 3 8 5 3 £ 3 8 0 9 D - 0 5 

2 

0.41 54989 1 28D 

01 

0.526532 05 09D- 05 

3 

0.41301 48775D 

01 

0.51 9 0 9 4931 6 D - 0 5 

4 

0.411521 633 OD 

01 

0.51 46566295D- 05 

5 

0.4107 03 1 778 D 

01 

0.51 23226637D- 05 

6 

0.41 092864 2 2D 

01 

.5120752653D-05 


ND= 1 FREQUENCY* 4 . 1 1 028D OORAD^SEC AMPLITUDE* 5.12006D-06 NIT* 6 


Figure 1 -18a . 
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LST SYSTEM-NUMERIRRL SOLUTION DF 1 +N^b 1+N^N<*G2=U . 

GAMMA= 1 .33 GOOD 07 TO* 1.0000 00-01 JV= 1.000000 05 -JG* £.100000 00 
KP= £.16 0 0 00 0£ KI= 9. 70 00 00 03 H* 6.0 000 00 02 K0= 5.758350 03 
Kl= 1.371020 03 K£= 1.000000-03 K= £.000000-01 


ITERATION 

FREQUENCY 

«ji 0 

0.411 028451 40 

1 

0.40320198810 

£ 

0.39804746390 

3 

0.39488 978 050 

4 

0.39295255330 

5 

0.39202408150 

6 

0.39226 047 060 


N0= 1 FREQUENCY* 3.92377 


AMPLITUDE 

01 0 .51 £00572570- 05 

01 0.494 09485700-05 

01 0.47991383650-05 

01 0.471 04448530- 05 

01 0.46567959110-05 

01 0.46314186740-05 

01 0 . 46286 028630- 05 

0 ORRL/SEC RflPL I TUOE= 



4.687660—06 M I T = 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N<>G1+N*N*G£=0 . 

GAMMA* 1.380000 07 TO* 1.000000-01 • JV* 1.000000 05 JG= 2.100000 00 

KP= £.160 0 00 02 * K 1= 9.7 0 00 00 03 H* 6. 00 GOOD 0£ K0= 5.758350 03 

Kl* 1.371020 03 K£= 1. 00000-03 K3= 3.000000-01 

ITERATION FREQUENCY AMPLITUD 


0 

0.39237667770 

01 

0.46276626930-05 

1 

0 .38386892930 

01 

0.44309736190-05 

2 

0.8 r 7 0977 f 0 1 D 

01 

0 .4271 9358320- 05 

3 

0.37318253090 

01 

0.41702164810-05 

4 

0.37067333130 

01 

0.41 0696904 00-05 

5 

0.37020428470 

01 

0.4 09 04 088770— 05 

6 

0.37 01 950 0790 

01 

0 .40847537970-05 

r 

0.37045117770 

01 

0.40846149610-05 


ND= 1 FREQUENCY* 3.704950 OORAD/SEC AMPLITUDE* 4.085730-06 NIT* 7 


LST SYSTEM-NUMERICAL S 
GAMMA* 1 .380000 07 TO* 
K P= £.160000 02 K I = 
K 1= 1.371020 03 K£* 


OLUTION DF 1 +N*G1+N*N*S£= 0 . 
1.000000-01 J V= 1 . 0 0 0 0 0 0 05 

9.7000 0 0 0 3 H * 6 . 0 0 0 0 0 0 02 

1 . 0 0 0 0 00“ 03 K3= 4.00 0 0 OD- 0 1 


Jb= £.100000 00 
K 0 = 5.75 8 3 5 0 0 3 


ITERATION 

0 

0 . 

FREQUENCY 
37 0494,8 1 260 

01 

AMPLITUDE 
0 . 4 0357268240- 05 

1 

0 . 

35906003020 

01 

0.38523260050-05 

2 

0 . 

35085157780 

01 

0.36572065410-05 

3 

0 . 

34516064410 

01 

0 . 35238 05 1 4 50-05 

4 

0 . 

34124206450 

01 

0 .84 8 6 7 5 6 3 4 8 0 — U 5 

5 

0 . 

34 164 9 35 4 3D 

01 

0 . 3 4 3 3 1691 8 5 0 - 0 5 

6 

0 . 

34 083290080 

01 

0.34187469390-05 

7 

0 . 

34122344610 

01 

0.34*84368570-05 

3 

0 . 

34 1 22943870 

01 

0.34171292720-05 

1 FREQUENCY* 3.41234D 

OORAr 

/SEC AMPL I TUBE* 3.41 64 1 0- 06 


Figure 1-1 8b. 



43 


LST SYSTEM-NUMER I CRL SOLUTION OF 1+N*G1+N*N^G2=0. 

GAMMA* 1 .38 0 0 0D 07 Td= 1.00 GOOD- 01 -JV= l.OOOOOD 05 Jb= 

KP= 5.160 0 OD 05 KI = 9.700000 03 H= 6. 00 GOOD 05 K0= 

Kl= 1.371020 03 K2= 1.000000-03 K3= 5.000000-01 

RMRLITUDE 
0.34164055030-05 
0 .3 0335459870—05 
0.5 7389S040 ID- 05 
0.551 079556 00-95 
0.53945954140-05 
0.57467997540-05 
0.31 533999 1 30- 05 
0. 56932383770- 05 
0.55470343270-05 
0.55458166550-05 
0.57504338340-05 
0 . 1 99 0339 0670— 05 
0.54089859510-05 
0.31 3 0 1 395480- 05 
0 . 3 0692755850- 05 
0 .26793983350— 05 
0.53245389560-05 
0.33835915440-05 
0 . 3 0576 061520— 05 
0 . 27656*3468 00— 05 
0.83949840450-05 

ND= 1 FREQUENCY* 2.941550 OORAD/SEC AMPLITUDE* 2.39h980-06 


LST SYSTEM-NUMER I CRL SOLUTION OF 1 +N+G1 +N-*N^G2= 0 . 

GAMMA* 1 .380000 07 TD= 1.000000-01 JV= l.OOOOOD 05 JG= 
KP= 8.160000 05 K I = 9.700000 03 H= 6.000000 02 K0= 
K 1 = 1.371 020 03 K2= 1.0000 OD- 03 . K3= 6.0000 00- 0 1 


ITERATION 

0 

FREQUENCY 

0.29418203530 

01 

AMPLITUDE 
0 . 5394984 0450- 05 

1 

0 . 3 1 688 1 42680 

01 

0.365719591 GD- 05 

O 

(m 

0 .870330761 ID 

01 

0 . 8 05299 0628 D- 05 

3 

0 .31605070450 

01 

0 . 2 8 6 3 9 0 9 2 8 9 D - 0 5 

4 

0 .26951249440 

01 

0 .2 01 3 0432350-05 

5 

0 .31 1554 4 8340 

01 

0 .27636719450—05 

t* 

0 .5633 8 u 9 7 3 6 D 

01 

0 . 1 937538 08 1 D- 05 

7 

U . 8 9 99865 32 1 0 

01 

0 .85679907640-05 

o 

0 . 19733495540 

01 

0.64 8 3 543 0 9 4 0 - 0 6 

9 

0.50535601840 

01 

0 9 0 9 7 6 6 3 0 480-06 

10 

0.855681 13600 

01 

0 . 15895889990-05 

1 1 

0.24850814170 

01 

0.15951056760-05 

12 

0 .29521825370 

01 

0.81 998709990— 05 

13 

0 .576 3 6 7 8 7 6 9 0 

01 

-0.131 92 1 65590- 04 

-5 .7636 

SO 00 -1.319 

580- 

05 

14 

15 

M . 9 9 9 9 9 9 9 HShfl 

0 0 

0 .99999^999611 00 

0.10 0 0 0 0 0 0 0 00 

01 

0.1 0 0 0 0 0 0 0 0 00 0 1 


NO* 1 FREQUENCY* 1.000000 OOP HD. -'SEC RMRLITUDE* 1.000000 00 


ITERATION 

FREQUENCY 


0 

0.34 1234 133 ID 

01 

1 

0.32410648200 

01 

2 

0 . 3 0773937 470 

01 


0 ■ 58 049952970 

01 

4 

0 .89805900980 

01 

5 

0 .59478965570 

01 

6 

0 . 3 28 3 5 9 9 9 3 4 D 

01 

7 

0 . 3 0397343 1 40 

01 

O 

•0.2321 4372560 

01 

9 

0.30015586400 

01 

10 

0.30977004720 

01 

1 1 

0.57495064330 

01 

12 

0.29430331640 

01 

13 

0.32957239420 

01 

14 

0 .31596592220 

01 

15 

0.3 0226 11 9770 

01 

16 

0 .D85893 0557D 

01 

17 

0.4684155030 

01 

18 

0 . 332 0392 03 OD 

01 

19 

0.3141 3885720 

01 

20 

0.29412203530 

01 


Figure 1 - 18c . 


2.100000 00 
5.758350 03 


N 1.1= 20 


2.100000 00 
5.758350 03 


HIT- 15 
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LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G1+N*N*G2®0 . 

GAMMA® 1 .33 CIO OB 0? TO= 1 . 00 GOOD- 01 JV= 1 . 000 COD 05 JG= £. 1 00000 00 

KP= 2. It. 00 OH 0£ K I = 9. 70 GOOD 03 H= 6.00000D 02 K0= 5 . 75835 D 03 

Kl = 1 .371 OSD 03 K£= 1.000 GOD-03 K3= 1 .0000 OH- 03 

I TERAT I ON FREQUENCY 

0 0 . 170000 OOOOD 01 

1 0 . 171 9691 394 D 01 

2 0. 173381721 OD 01 

3 0. 174 165639 OH 01 

4 0. 174589 163 ID 01 

5 0 . 1743 06553 7D 01 

6 Q.174S141543D 01 

ND= 1 FREQUENCY® 1 .74S23D OGRAB^SEC AMPLITUDE® 4.45254D-Q7 NIT® 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G1+N*N*G2=0 . 

GAMMA® 1 .3300 OB 07 T0= 1.0000 OB- 01 JV= 1.0000 OB 05 

KP= £.1600 QD3 02 ' K I = 9.7 00 0 OD 03 H= 6 . 0 0 0 0 OB 02 

Kl = 1 .37 1 OSD 03 K2= 1.00000B-03 K3= 1.00000B-02 

4 

ITERATION FREQUENCY 3 AMPLITUDE 

0 0 . 1 74822728 OD 01 0 .4452536775H-06 

1 0. 1752126778D 01 0 .4502796250B-06 

£ OM 75365 0190B 01 0 .4529774 839H- 06 

3 0.1754421 16 ID 01 0 . 4544379 007D- 06 

4 0 . 1 7548888 5 9D 01 0 .4552 1 38844D- 06 

5 0. 17551 72 026D 01 0 .4556845266B- 06 

6 0.1 755362634D 01 0 .4559779957B- 06 

ND= 1 FREQUENCY® 1 .75550B OORAD^SEC AMPLITUDE® 4.56178D-07 NIT® 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N+G 1+N*N*G£=0 . 

GAMMA® 1.33 00 OD 07 TO® 1.0 OOOOD- 01 JV= 1.0 OOOOD 05 J6= 2.1 OOOOD 00 

KP® 2 . 160 OOD 02 KI= 9.7 OOOOD 03 H= 6. 000 ODD 02 K0= 5.75335D 03 

Kl= 1.371 02D 03 K£= 1.00000D-03 K3 = 1.0 00 OOD- 01 

I TERAT I DN FREQUENCY AMPL I TUBE 


0 

0.1 75549719 ID 

01 

0 . 4 56 1 78 4 994 D — u6 

1 

0.179 

1 779346D 

01 

0.51 0161 2925D- 06 

cl 

0.181 

37 1 221 5D 

01 

0 . 542 0 1 84 2 1 4 D- 06 

•-v 

o 

0.132608181 ID 

01 

0 . 559583 0487D- 06 

4 

0 . 1 83 

£62 172 ID 

01 

0 . 5 6 8 7 8 9 7 9 0 4 D - 0 6 

5 

0.182 

5348001 D 

01 

0.57 3 8 8 0 8 5 6 0 1* - 0 6 

6 

0.132 

642 032 OD 

01 

0 .576S742238D-06 

( 

0.182 

7338 002D 

01 

0.5780171 524D- 06 

y 

0 . 1 82 

7323626D 

01 

0 .57891 074 07D- 06 

9 

0.182 

3 17472 ID 

01 

0 .5793698723D-06 


ND= 1 FREQUENCY® 1 .83840D OORAB-'SEC AMPLITUDE® 5 . 79735D- 07 NIT® - 9 


JG= 2.1 OOOOD 00 
KG® 5 .7583 5 D 03 


AMPLITUDE 
0.45000000 0 OD- 06 
0.4377777686B-06 
0 .440 1 73 0766D- 06 
0 . 4 4 26 0 1 S3 22D — 06 
0 .444 1 635449D-06 
0 .444781 2824D- 06 
0.4451 726996D- 06 


Figure 1-1 9a. 
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LST SYSTEM-NUMERICAL S 
GAMMA* 1 .38 00 OB 07 TD= 
KR= £.1600 01' 0£ K I = 
■ K 1= 1.371 OSD 03 K£= 


□LUTIDN UF l+^61+h^G£=0. 

1 . 000 OOP- 01 JV= l.OOOOOD 05 

9.7000 OD 03 H= 6. 0000 on 02 

1.0000 OD- 03 K3= £.00 0 0 OD- 0 1 


•J6= £. 1 00000 0 0 
K0= 5-75835D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 18384 0 0 034D 

01 

0.579735429 ID- 06 

1 

0 . 1886569754D 

01 

0.651 17421 56D-06 

O 

u 

0 . 19174334 03D 

01 

0 .695722 0£ 57 D— 06 

3 

0. 1936546542D 

01 

0.721 1874 01 ID- 06 

4 

0.1 92^9454.8 OD 

01 

0.7 4994 1 8438B * 06 

5 

0 . 1 933822886 D 

01 

0 .746801 0454D-06 

h, 

0 . 1 9475 1 5 1 77D 

01 

0.75 00348 03 ID- 06 

*** 

r 

0.1 95£344£39D 

01 

0.7521 033896D- 06 

8 

0.1 954 077725D 

01 

0. 7522564 047D- 06 


MD*. 1 FREQUENCY* 1 .95414D OORAD^SEC AMPLITUDE* 7.5£93£D-07 NIT* 8 


LST SYSTEM-NUMERICAL S 
GAMMA* 1 .380 ODD 07 TO* 
KF- £.160000 0£ KI= 
Kl= 1.37 10£D 03 K£= 


□LUTIDF CF 1+N*G1+N*N*6£=0 . 

1 .0000 OD- 0 1 JV= 1 .0000 on 05 

9 . 7 0 0 0 on 03 H* 6.0 0 0 0 OD 02 

1 . o 0 0 0 OD- 03 K3= 3 . 0 0 0 0 OD- 0 1 


J6= 2.10000D 00 
K 0= 5 . 75335D 03 


ITERATION 
3 0 
1 


4 

5 

6 
7 


FREQUENCY 
0 . 1954136££5D 01 
0 . £ 0 1 5669£59D 01 
0 . £ 0598 1 2684D 01 
0 . £ G638££449D 0 1 
0 . £ 0845 7323 5D 01 
0 .£0964141 lull 01 
0.21 0254214 ID 01 
0 .£1 04071397D 01 


AMP I TUBE 
0 . 7529324 744 D- 06 
0 . 843724989 OD— 06 
0 . 9 0342 1 £ 1 99D- 06 
0 . 962 08 07 4 53 D— 06 
0 . 971 922544 OD- 06 
0 . 9788£3 0989D- 06 
0 . 98228 0076 ID- 06 
0 . 98239542 1 9D- 06 


ND= 1 FREQUENCY* 2.10420D OORAD^QEC AMPLITUDE* 9.83224B-05 NIT* 7 


Figure 1-1 9b . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^61+N*N*S2=0 . 

GAMMA* 1 .3800CD 07 TO* 1 . 0 000 00-01 JV= 1.00000B 05 JG* 2 . 10 00 OB 00 
K P* 2. 16 00 OH 02 KI = 9.70000P 03 H= 6. 00 00 OB 02 KO* 5.75S35D 03 
K 1 = 1.37 102H 03 K2= 1 . 0000 CD- 01 K3= 1 . 000000-01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .40 0000000 OB 

01 

0.50 0 0 0 0 0 0 0 0 B - 0 5 

1 

fj . 4 06'9u4y4 y 9B 

01 

0.5 06 08167590- 05 

£ 

0.4 09 6 U U2875D 

01 

0 . 5 0 75 1 4 697 OB— 0 5 

3 

0.41 053525 05D 

01 

0. 5 0770691 71 B- 05 

4 

0.4 10994741 OB 

01 

0 ■ •_* 0 1‘ b.o Jc- 934 B— iJ-_i 


ND* 1 FREQUENCY* 4.10996B OQRAD^SEC AMPLITUDE* 5. 0763 OB- 06 NIT* 4 


LST S 
GAMMA 
KP= 2 
Kl= 1 


YSTEM-NUMERI CAL S 
= 1 .33 00 OB . 07 TO* 
.16 0 0 0 D 02 K I = 

.37102B 03 \ K2= 


□LUTIDN DF 1+N^G1+N*N*G2=U. 

1 . 0000 OB- 0 1 J'v - 1 . 0 0 0 0 OB 05 

9 .700 0 0 B 03 H= 6 . 0 0 0 0 0 B 02 

1 . OOuOOB— 01 K3= 2 . 0 00 00B-01 


JG* 2. 1000 OB 00 
KO* 5.75835D 03 


I TERAT I DM FREQUENCY 

0 0.41 09963752B 01 

1 0 .403 17097 1 2B 01 


2 0.398 00 037 04B 01 

3 0. 39482375 13B 01 

4 0. 392865 1596B 01 

5 0.3 9 2 0 0 0 0 3 7 0 B 0 1 

6 0 .392229276r« B 01 


AMPLITUDE 
0 .5 076302789B- 05 
0. 4897281 315D-05 
0 .475535698 1 D- 05 
0 .46664 3 861 1 B — 0 5 
0.461 252 0582B- 05 
0 . 45385 1 2642B- 05 

0.4585636 143B-05 


NB= 1 FREQUENCY* 3.92343B OORAD^SEC AMPLITUDE* 4.58466B-06 NIT* 6 


LST SYSTEM-NUMERICAL SOLUTIDN DF 1+N^G1+N^N*G2=0 . 
GAMMA* 1 .38 00 OB 07 T0= 1.00 00 OB- 01 JV= 1 .000 0 OB 05 

KP= 2. 16 GOOD 02 KI* 9.70000B 03 H= 6. 0 0 00 OB 02 

Kl= 1.37102B 03 K£= 1.00000D-01 K3= 3.00000B-01 


JG* 2.10000B 00 
K0= 5.75835D 03 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .3923434 00 6D 

01 

ll .458466259 1 B- 05 

1 

0.383 0573 188D 

01 

0 . 4 3 8 3 9 8 6 1 0 5 B - 0 5 

2 

0.3767396 17 OB 

01 

0.422171 1154D-05 

3 

0 .3727 07392 ID 

01 

0.41 174 951 78B- 05 

4 

0.3701 052335B 

01 

0 .4052404 157B- 05 

5 

0 . 3 697 3 3 4 2 9 3 B 

01 

0.4 0373892 04D- 05 

6 

0 . 369 6 4 5 1 o 4 5B 

01 

0.4 03 0788857B- 05 

r 

0 . 37 0 0278587B 

01 

0.4 03 179332 OB- 05 

8 

0 .370 068438 OB 

01 

0.4 034395 196B- 05 


ND= 1 FREQUENCY* 3.70157D OORAB/SEC AMPLITUDE* 4. 0361 ID-06 NIT*. 8 


Figure l-20a. 
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REPRODUCIBILITY of the 
ORIGINAL PAGE is POOR 


LST SYSTEM-NUMER I CfiL SOLUTION OF 1+N*61+N*N^G£*0 . 

GfiMMfi* 1 . 33 0 0 OB 07 TO= 1 .0000 OB- 0 1 JV= 1 . 0 0 0 0 OB 05 -JG= 2.1000 OB 0 
KP= £.16 0 0 OB 02 KI= 9 .7 000 OB 03 H* 6. 0000 OB 02 K0= 5.75835B 0 
Kl= 1.37102B 03 K£* 1.00 00 OB- 01 K3= 4. 00 GOOD- 01 


ITERATION 

FREQUENCY 



fiMPLITUBE 

0 

0 . 37 01567 1 71 B 

01 

0 .4036108745B-05 

1 

0. 35828 J6527D 

01 

0. 

3793426426B- U5 

iu! 

0 . 3 4965 5 3 1 6 8 B 

01 

0. 

358954 5285B- 05 


0. 34355 U973B 

01 

0 . 

344778 07 17B- 05 

4 

0 .339293 003 SB 

01 

0 . 

3 354372236B- 05 

5 

0 . 3398 1 93 057B 

01 

0. 

335229799611-05 

6 

0 .3392 1 47326B 

01 

0 . 

334 06 3 1 4 74 L - 0 5 

7 

0 .33951 77657B 

01 

0. 

333969 1379B- 05 


NB= 1 FREQUENCY* 3.39539B OORAD^SEC fiMPLITUBE* 3.33S65B-06 .NIT* 7 


LST SYSTEM-NUMERICAL SOLUTION OF’ 1 +N*G 1 +N*N^G£= 0 . 

GAMMA* 1 .33 00 OB 07 TO* t.OOOOOB-Ol JV= 1.00000B 05 Jb= £.1000 OB 00 
KP = £ . 1 6 0 0 OB 02 K I * 9.70 0 0 OB 03 H* 6.00 0 0 OB 02 K 0= 5 . 75835B 03 
Kl = 1.37102B 03 K2* 1.00000B-01 K3= 5.00000B-01 


ITERATION FREQUENCY 

0 0 .33953S5263B 

1 0.32091 1352 OB 

2 0 .3003614935 B 

3 0.1 1354 1332 IB 

1 1.18541B 00 -8.49 


fiMPLITUBE 

01 0 . 333865 4 6 74 B - 05 

01 0 . 29762 1 497 6D— U5 

0 1 0.2 5 6 0 3 1 3 8 6 5 B - 0 5 

01 - 0 . 84973759 04B- 06 

38B-Q7 


4 0.100 0 0 0 0 098 B 01 0 . 9999999 03SB 0 0 

5 0 . 1 0000 00 00 OB 01 0 . 1 OOOOOOOOOB 01 


ND= 1 FREQUENCY* 1 . 00000B OORfiB-'SEC fiMPLITUBE* 1.U0000B 00 NIT* 5 


Figure 1 -20b . 


CO o 
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LST SYSTEM-NIJMERICfiL SOLUTION OF 1+N^G1+M*M*G2=0 . 

GAMMA* 1 .33 00 OH 07 10= 1.00000P-01 JV= 1.000 OOP 05 J6= 2. 10 00 CD 00 

K P = £.160 0 0 D 0 2 K I = 9.70 0 0 0 D 0 3 H = 6.0 0 0 0 0 D 0 £ K 0 = 5 . 7 5 8 3 5 D 0 3 

K 1 = 1.371 02B 03 K£= 1 .00 0 0 OP- 0 1 K3= 1 . 0 0 0 0 OP- 0 1 

ITERATION FREQUENCY AMPLITUDE 


0 

0.1 

3 0 0 0 0 0 0 0 OP 

01 

0 . 6 0 0 0 0 0 0 0 0 OP- 06 

1 

0.1 

329985273D 

01 

0 . 5796322434 D- 06 

£ 

0.1 

35465 0 6 8 3 P 

01 

u .587191 6572P- 06 

3 

0.1 

3683794 08P 

01 

0 .59 3 0 4 8 6 3 5 4 P — 0 6 

4 

0.1 

375561 £88D 

01 

0 . 5 9 6 3 £4074 8 D - 0 6 

5 

0.1 

379253743D 

01 

0 .59 78864 735 D- 06 

6 

0.1 

379958 0£7B 

01 

0 .59786 0951 2D- 06 


NB= 1 FREQUENCY* 1 .83003D OORAD/SEC AMPLITUDE* 5.98399D-07 HIT= 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G1 +N*N*G£= 0 . 

GAMMA* 1 .330 OOP 07 TO= 1.00000P-01 JV= 1.00000D 05 JG= S.10G00D 00 
KP= £.16 00 OP 02 KI = 9. 70 00 OP 03 H= 6. 000 OOP 02 K0= 5.75835P 03 
K 1 = 1.37102P 03 : i K£= 1 .00 00 OP- 01 K3= £.000 OOP- 01 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 . 188007791 7D 

01 

0 .5983992937D- 06 

1 

0.1 9294 1 2767P 

01 

0. 672325 1343D- 06 

U 

0 . 196 1241 662P 

01 

0.71 86786352P- 06 

3 

0 . 1981 174432P 

01 

0.7452828246D-06 

4 

0.1 979638 157D 

01 

0.771 5 3345 7 9 B— 06 

cr 

0 . 199 007872 OP 

01 

0 . 7 7 4 7 9 3 7 8 6 4 P — 0 6 

6* 

0 . 1996 13 0859 P 

01 

0 . 777426886 1 P-06 

“5 

1 

n . 1999 099 796 p 

01 

0. 7783768 185P- 06 

8 

0.1 9992 07 66 1 P 

01 

0.7786887285D-06 


NP= 1 FREQUENCY* 1 .39932P OORAP'SEC AMPLITUDE* 7.7S826P-07 HIT= 8 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N+G 1 +N*N^G2= 0 . 

GAMMA* 1.38 00 OP 07 TO* 1.00 00 OP- 01 JV= 1.00 00 OP 05 JG= £.1000 OP 00 

KF‘= 2.1 6 0 0 0 P 0 £ K I = 9.7 0 0 0 0 D 03 H = 6 . 0 0 0 0 0 P 02 K 0 = 5 . 7 5 8 3 5 P 0 3 

K 1 * 1 .371 02D 03 K£* 1 . 0 0 0 0 OP- 0 1 K3= 3 . 0 0 0 0 OP- 0 1 

ITERATION FREQUENCY AMPLITUDE 


0 

0.1 

99932 3 8 76P 

01 

0 . r »■ 8 8 262524 P U 6’ 

1 

0.2 

062672759D 

01 

0.87343191 04P- 06 

o 

L. 

0.2 

109122 09 1 P 

01 

0 .93622496 01 P-06 

* 

0 .2 

1 15867393P 

01 

0.10 0792 0 039D- 05 

4 

0 . 2 

13285322 ID 

01 

0 . 1 0 0 9 0 06 2 04 D - 05 

5 

0 . 2 

14 566567 6 P 

01 

0.1016 0 09 078P- 05 

6 

0 .2 152649862 D 

01 

0 . 1 02 0 05 1 049D- 05 

T 

t 

0.2 

1550 IS 07 ID 

01 

0.1 02 0575884 P- 05 

8 

0 . 2 

155123 139D 

01 

0.1 02 1 438689D-05 


ND* 1 FREQUENCY* 2.15534D OORAD'SEC AMPLITUDE* 1.02156P-06 NIT* 8 


Figure 1 -21a 
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LST S V STEM- MU M E R I C RL SOLUTION OF 1+N^G1 +N*N^G2®0 . 

GHMMA = 1 . 38 0 0 OIi 07 TO= 1.0000 0 D- 0 1 J V= 1 . 0000 OD 05 .JG= £.1 0 0 0 OD 0 0 

KP= £.160000 0£ K I = 9.700000 03 H= 6.000000 0£ K 0 = 5.753350 03 

K 1 = 1.371 0£D 03 K8= 1.0000 00 - 0 1 K 3= 4.0000 OD- 0 1 

AMPLITUDE 
0 • 1 0£'l 5637470- U5 
0.1 1539385170-05 
0.1 £6 137 0534 0-05 
0.1 ££574536 ID- 05 
0.1 £9 09363 190-05 
0 . 1 343 039 1 530- 05 
0.13737431730-05 
0.1 33869 11 £00- 05 
0 . 1 393 1 66 0 030- 05 
0 . 1 3 9 3 8 9 3 7 9 5 0 - 0 5 
0 .13947084940-05 

N0= 1 FREQUENCY® £.334350 QGRAD/SEC AMPLITUDE® 1.395670-06 NIT® 10 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G1+N^N*6£® 0 . 

GRMMA® 1 .380000 "07 TO 1 .0 000 00-01 JV= 1 .0 000 OD 05 

KP= £.160000 0£ KI® 9.700000 03 H= 6.000000 03 

Kl = 1.371030 03 K3= 1.000000-01 K3= 5.000000-01 

ITERATION FREQUENCY AMPLITUDE 

0 0.33843540030 01 0 . 1 39566676 OD- 05 

1 0 .35356 09 9 8 8H 01 0 . 1637 1938 170-05 

3 0 . 363 0 1 3974 1 D 0 1 0.1 9 0988457 1 D- 05 

3 0 . 3 3 063669 OED 0 1 0 . 3 1 798339840- 05 

4 0 .39793 1371 8D 01 0 .346965431 3D- 05 

5 0 . 3937 0 1 05750 01 0.351 £6983340-05 

6 0.35393377510 01 0.17739539830-05 

7 U . 3693 3 657630 01 0.1 94 1 4483350— 05 

8 0 . £79 0768 0 070 0 1 0.31 056779390-05 

9 0 . 38 1 633 09360 0 1 0 . 38933476740- 05 

10 0.38141 449450 0 1 0 . 893556578 OO- 05 

11 0.30159941590 01 0.35750971960-05 

18 0.18197513040 01 0.33756674810-06 

1 3 0 . 1 844 1 0 06 140 0 1 0 . 5337 1 498730- 06 

14 0.1 9979978560 0 1 0 . 7845 1 7 08660- 06 

15 0.316 03553990 01 0.1 0368854 070- 05 

16 0 .83164374770 01 0 . 1 88 1 65893 00- 05 

17 0.34734571630 01 0.15801503140-05 

1 3 0 . 85 034933380 01 0.18 07034 1 3 OD- 05 

19 0 .3643365 1 1 80 0 1 0.1 9539 1 048 OD- 05 

20 0 .88190369760 01 0 .831 75365580-05 

HD® 1 FREQUENCY® £.819090 OORAD/SEC AMPLITUDE® 8.317530-06 HIT= £0 


JG® 8.100000 00 
K 0® 5.758 3 5 0 U 3 


ITERATION 

FREQUENCY 


0 

0.81553398610 

01 

1 

0.33448351 9 ID 

01 

8 

0 .80353703850 

01 

3 

0.31803474800 

01 

4 

0 . ££7 1 1386140 

01 

5 

0 . £8£8 , ^ c isx , -H OD 

01 

6 

0 .£361 0756840 

01 

7 

0.33771 843 15Ii 

01 

3 

0 .33384993330 

01 

9 

0 .83330506930 

01 

10 

0.23337717140 

01 


Figure 1 -21 b 
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LST SYSTEM-NUMERICAL SOLUTION QF 1+N^G1 +N^N<*G2*0 . 

JAMMA® \ .IllodlOfi 0? T0« \ , 000005- Pi JV® 1 . g 0 0 u 00 05 JS* 2. 10 00 OIi 00 

kp=» a.isooori oa ki= 9.700000 03 h= s.oooood 02 ko= 5.753350 03 

Kl= 1.37 1020 03 K2= 1.00 GOOD 00 K3= 1.000000-03 

HERAT I DM FREQUENCY 

0 0.430 00 0 0 0 0011 01 

1 0.43 0523 1336D 01 

2 0.43 0528335 OD 01 

3 0.43 0874487 OD 01 

4 0.43Q8634353D 01 

ND= 1 FREQUENCY* 4.3G865D OORAD-'SEC AMPLITUDE* 5.24185D-06 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N^G2=0 . 

GAMMA* 1 .330 ODD 07 TO* 1.00 GOOD- 01 JV= 1.00 GOOD 05 JG= 2 . 1 0000D 00 

KP= 2.160000 02 K I = 9.7000 OD 03 H= 6. 00 GOOD 02 K0= 5.75835D 03 

Kl= 1 .371 02D 03 4 K 2= 1.00000D 00 K3= 1.00000D-02 

AMPLITUDE 
0.5241 85 0229D- 05 
0.522 0253 01 3D- 05 
0.5 1944726 04H-05 
0.5194433793D-05 
0.51 9386774 9D- 05 

ND= 1 FREQUENCY* 4.29161D OORAD-'SEC AMPLITUDE* 5.19344D-06 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION DF 1+N*G1+N*N*G2=0 . 

GAMMA* 1.33 GOOD 07 TO* 1.00 GOOD- 01 -JV= 1.0000GB 05 

KP= 2 . 1 6 0 0 OD 02 K I = 9.7 0 0 0 OD 03 H* 6.00 0 0 OD 02 

Kl= 1.37102H 03 K2= 1.0 000 OD 00 K3= 1.0 00 OOP- 01 

I TERAT I DM FREQUENCY AMPL I TUPE 

0 0.429 16 0758 OD 01 0 .5 1 9344 0699D-05 

1 0 . 4222646483D 01 0.5 0262 Q0492D- 05 

2 0.4 177 09477 OD 01 0 .4895733868D-05 

3 0 .4 149095 0040 01 0 . 48 1 4 07 1 54 OD- 05 

4 0.41-31 739544 D 0 1 0.47 646 U78 r 1 OD— U5 

5 0 .4 1 262777 7SD 01 0 . 4746822755D-05 

6 0.41281 02 003D 01 0 . 4743764 1 69D- 05 

ND* 1 FREQUENCY* 4.12899D OORAD-'SEC AMPLITUDE* 4.74269D-06 NIT* 6 


Jb* 2.1 000 OD 00 
K0= 5 . 75 835 D 03 


ITERATION FREQUENCY 

0 0.43 08645 193D 01 

1 0.429859 05 14D 01 

.2 0 .423729496-4P 01 

3 0 .42919135 05D 01 

4 0.4291 7624 17D 01 


AMPLITUDE 
0.52000 0 0 0 0 OD- 05 
0 .52283 039 0 1 D- 05 
0. 52224534390- 05 
0 .524 5748 72 8D- 05 
0.5243241 015D-05 


Figure 1 -22a . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^61 +N^N*62=0 . 


GAMMA 3 1.38 00 OB 

07 TO 3 

1 

. o 0 0 0 0D- 

01 

JV 3 

1 . 00 GOOD 05 

.JG- 2.10000B 

00 

KP= 2. 16 GOOD 02 

. K I = 

9 

.70 0 OOD 

03 

H= 

6.00 0 0 OD 02 

KO 3 5.75335D 

03 

Kl= 1.37102D 03 

K2= 

1 

. 00 GOOD 

00 

K3= 

2. 000 ODD -01 




ITERATION FREQUENCY 

0 0 .4128986601 D 01 

1 0.4 0 3 8 3 3 3 3 1 8 D 0 1 

2 0 .3377 0780421) 01 

3 0 .3933 1 4 08 9 90 01 

4 0.39 13 0588350 01 

5 0.33 0341 7445 B 01 

6 0 • 39 U3686865D 01 

7 0.331 0473948D 01 

ND= 1 FREQUENCY 3 3.91 085D OORAD'SEC AMPLITUDE 3 4. 175590*06 NIT= 7 


AMPLITUDE 
0.4742691 1 03D-05 
0 .45344537 3 3 B — U 5 
0 .43699 1 6734B— 05 
0.4264565926D-05 
0.41 93377967B- 05 
0.4183451 03611-05 
0 . 4 1768 17051 B— 05 
0.4 17557537 OB- 05 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N^G 1 +N^N*G£ 3 0 . 

GAMMA 3 1 .38 00 OH ^07 TO= 1 .00000D-01 -JV= 1 . 00000D 05 JG= 2.10000D 00 
KP= 2. 16 0 0 OD 02 K I = 9.700000 03 H= 6.00000D 02 ' K0= 5.75835D 03 
Kl= 1.37102D 03 K2= 1.00000D 00 K3= 3.00000D-01 


ITERATION FREQUENCY 

0 0 . 39 1 085 0693 D 0 1 

1 0 .3796 023 096 D 01 

2 0 . 37 1 28 03 035 D 01 

3 0.36541 73534 D 01 

4 0 .361 39855 77D 01 

5 0.361 942 0497D 01 

6 0.3611 493 035D 01 

7 0 .361 5626 743 B 01 

8 0.361 56 19356D 01 


AMPLITUDE 
0.41 75586358D- 05 
0 . 392 6 7344 4 4 B- 05 
0.372 054 1479D- 05 
0.357766594 ID- 05 
0 . 3485 1 8 7376 D- 05 
0 . 3483 1 88329D- 05 
0 .3468979466B-05 
0 .346872725 o I) — U 5 
0. 346741 56 19D- 05 


ND= 1 FREQUENCY 3 3.61565D OORAB/SEC AMPLITUDE 3 3.46667B-06 NIT= 3 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +N*N>G2 3 0 . 

GAMMA 3 1.38 GOOD 07 TO 3 1.00000D-01 JV 3 1.00 00 OB 05 JG 

K.p= 2. 16 00 OB 02 KI = 9. 7 0 GOOD 03 H= 6.0 000 OD 02 KO 

Kl= 1 .371 02 D 03 K2= 1.00 GOOD 00 K3= 4. 00 GOOD- 01 

ITERATION FREQUENCY AMPLITUDE 


0 

0 .36 15649 056B 

01 

0 

. o 4 666 78 496 D — U5 

1 

0 .343429794 9B 

01 

0 

.3 093 039 01 2D- 05 

L. 

0 .324 067764 9D 

01 

0 

. 2 6 69 93 0 1 83D - 0 5 


3 0. 13 09721 258D 01 -0 . 1 572573 193B- 06 

1 1 .30972D 00 - 1 .572570-07 

IHC900I EXECUTION TERMINATING DUE TO ERROR COUNT FOP ERROR NUMBER 261 


3 2 . 1 OuOuD 00 
3 5.7 5 8 3 5 B 0 3 


Figure 1 -22b . 
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LST SYSTEM-NUMFR ICAL SOLUTION OF 
GAMMA 3 1.380C00 C7 l.OOOOOD- 

1+N*G1+N*N*G2=0. 

01 J V= l.OOOOOD C 5 

JG 3 2.10000D 

00 

K P= 2 . 160 OOD 02 KI= 9.70000D 

Kl 3 1.37102D 03 K2= l.OOOOOD 

03 H 3 6 • OOOOOD 02 

00 K 3= 1 • 00000 D- 03 

K 0= 5.75835D 

03 

IT EP AT I ON FREQUENCY 

0 0.21000000000 01 

1 0.21058192930 01 

AMPLI TUOE 
0.5800000000D-06 
0.576440 8218D-06 



2 0.21 10683445D QL 

3 0.21 13400637D 01 

4 0.2113458576D 01 

0.57779394370-06 
0.5784885398D-06 
0.578755 19580-06 




ND 3 1 FREQUENCY 3 2*113520 QORAD/SEC AMPLITUDE 3 5*788290-07 NIT- A 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1+N*N*G2=Q. 

GAHMA= 1.38QfQO C 7 TO 3 1,000000-01 J V= 1*000000 05 J G= 2,10000D 00 

K P- 2,160000 02 K I- 9.70000D 03 H- 6,000000 02 K0= 5.75835D 03 

K 1= 1.371020 03 \ K2= l.OOOOOD 00 K3 3 1.00000D-02 


ITERATION 

. FREQUENCY 


AMPLITUDE 

0 

0.21 13522368D 

01 

0.5788294272D-06 

1 

0. 21 180918760 

01 

0.5 850930129D-06 

2 

0.21 1989 5 2 1 8D 

01 

C. 58873250600-06 

3 

0.21 2090 5 1 68D 

01 

0.5907279733D-06 

4 

0, 21 21 5 52 4440 

01 

. C. 5S18624186D-06 

5 

0.21 219679300 

01 

0*552602 31040-06 

6 

0.21222618840 

01 

0. 5930924560 D- 06 


ND= 1 FREQUENCY 3 2.12248D OOR.AD/SEC AMPLITUDE 3 5, 93446D-07 NIT= 6 


LST SYSTEM-NUMERICAL SOLUTION OF 1 + N*G U N*N*G2 3 0 . 

GAMMA= 1.380000 07 TO 3 1. 000000-01 J V= l.OOOOOD 05 JG 3 2.100000 00 

KP= 2*160000 C2 K I = 9.70000D 03 H= 6,000000 0 2 K 0= 5.75835D 03 

K 1= 1.371020 03 K2= l.OOOOOD 00 K3= l.OQOOOD-Oi 


ITERATION 

FREQUENCY 


AMPLI TUD: 



0 

0.2122476046D 

01 

0.5 93446 0004D-06 



1 

0.2167383373 D 

01 

0.6662 85 8378 D- 06 



2 

0. 21 95828C91D 

01 

C. 711400 74880-06 



3 

0.22 1 283 2 34 2D 

01 

0.7 3721802100-06 



4 

0.2224826790D 

01 

0.74964577550-06 



5 

0.2227666826D 

01 

0.759949 1037D-06 



6 

0.22301331990 

01 

0.786206 7825D-06 



7 

0.22267340260 

01 

0.7706783119D-06 



8 

0.22303627910 

01 

C. 771421 38800-06 



9 

0.223 1377553D 

01 

C *77085962860-06 



ND= I FREQ 1 

JENCY = 2 • 2 3 1 4 ID 

OQRAO/SEC AMPLITUDE 3 7.71332D-07 

NIT- 9 



Figure l-23a. 
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1ST SYSTEM-NUM C R ICAl SOLUTION OF l + N*Gl + N*N*G2-0. • _ 

G AMM A* 1.380000 07 TO* i.OQOOOD-OL JV= 1*000-000 05 JG= 2.100000.00 

KP = 2.160000 02 KI* 9.700000 03 H* 6.00000D 02 K0= 5.758350 03 


Kl= 1.371020 


K2 = 1 .000 000 


K 3 = 2.000000-01 


ITERATION FREQUENCY AMPLITUDE 

0 0.223140651 90 01 0.7713 32 3 18 2Q-Q6 

0.22918969950 01 C . 8 703 78 97 240-06 


2 

0.23332686360 

01 

C .93602606470-06 

3 

0.23674924000 

01 

C. 97474662960-06 

4 

0.23 75837684D 

01 

0.1C0150 89090-05 

5 

0.23797925870 

01 

0. 10134590340-05 

6 

0.23822084950 

01 

0.1C2166187BD-05 

7 

0.23026721340 

01 

0.1C 26 021014 D-05 

8 

0.2384497 1350 

01 

0.1C26574826D-05 

9 

0. 23847699590 

01 

0. 1C279821240-05 

10 

0.23852289 8 7D 

01 

0.10283023410-05 


ND= i F REOU ENCY= 2.385480 OORAD/SEC A MPLI TUDE = 1. 028880-06 NIT= 10 


j L ST SYSTEM-NUMERICAL SOLUTION OF 1 + N*G1 + N*N* G2=0. 

GA MMA = 1 .38000D 07 TO* 1.00000D-01 JV= l.QOQO OD 05 JG= 2.10000D QQ 

' KP= 2.16000D C2 KI= 9.7C000D 03 H = 6.00000D 02 K0= 5.758350 03 

Kl= 1.371020 03 K 2= l.OOOOOD 00 K3= 3.000000-01 


ITERATION FREQUENCY AMPLITUDE 

0 0.23 8598396 10 01 0. 1 C28 88 15 16D-05 

1 0. 296780 i 350 0 01 0.11650979630-05 

2 0.25356528850 01 0 • 1 266 39 8 36 2D-0 5 



0.255 78950900 01 0. 1 2 92 37 439 80-05 
0.257 76 36773D 01 0 . 1 39420 10290-05 
0.2588264 309D 01 0. 136989 89180-05 
0.25948114380 01 0 . 1 3 84 32 1 5460- 05 
0 . 2599548 808D Ql 0.139525 018 10-05 
0.26029329430 01 0.14027589710-05 
0.2605653351D 01 0. 14 10 3999 390- 05 
0.2605424268D 01 0 . 1 425 49 9 8 1 8 D-05 
Q. 26 12458009D 01 0. 1 4290249370-05 
0.26151963730 01 0 . 1 4 29 3 5 34 9 80-05 
0.26 152382Q7D 01 0 . 1 430 13 397 8 D-05 


NO* 


FREQUENCY* 2.615340 OORAD/SEC AMPLITUDE* 1.43025D-06 NIT* 13 



Figure 1 -23b . 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G 1+N*N* G2 = 0. 

Qb MMA* 1.3B00QQ Q7 T 0 = l.OOOOQD-Ql J V= l. QQOUQD 05 JG= 2.100000 OC 

KP= 2.160000 02 KI* 9. 70000D 03 H= 6.000000 02 K0= 5 . 758350 02 

Kl= 1.371020 03 K2= 1. 000000 UO K 3= 4.000000-01 

iteration FREQUENCY AMPLI TUDt : 

0 0.261 533 8.32 90 01 0.143025 222 70-05 

1 0.2749 1451770 01 0.16 57689676 D-05 

2 0. 271550C 806D 01 0. 1 6331919630—05 

3 0.28568670450 01 0. 1 992 147664D-05 

4 0.30181244280 01 0 . 2 2 43 11 3 599P-Q5 

5 0.32362731020 01 0. 2 63005 92 730-05 — 

6 0.2985502 1500 01 0 . 2 252 5 79 94 10-05 

7 0.317688 3 1 180 01 0.256 75 3901 2.0- 0 5 

8 0. 19 141 1 23C9D 01 0. 4 2841 346750-07 “ ” 

9 0.20 8338 6846D 01 0.6 4 30 28 84600-06 

10 Q.22 882432Q30 01 0.935 746 10Q8D-06 

11 0. 24575 10229D 01 0.11954391140-05 

12 0. 25'99 1 5003 2D 01 0.14282643 080-05 

13 0. 271976961 70 01 0* 1635 828728 0-0 5 

14 0* 28 2440 43 39 D 01 0 . 1 E2 099 973 8 D-05 ~ 

15 0.29220335750 01 0. 1 99 380 0 62 30-05 

16 0.2P 804 18289D 01 0 . 2 C 7 1 62 534 9D-Q5 

17 0. 30 1537490C0 01 0 . 2 2590 8 07 78 D-05 


18 

0.32966725560 

01 

0.27652305310-05 

19 

0. 30 0 5244 2540 

01 

0.22179562590-05 

20 

0. 32077926G7D 

01 

0. 25743091210-05 


ND= 1 FREQUENCY* 3.20779D OORAD/SEC AMPLITUDE* 2.574310-06 NIT* 20 


LST SYSTEM-NUMERICAL SOLUTION OF 1+N*G1 +N*N* G2=0 . 

GAMMA* 1.33CCC0 Q7 TO* 1.000000-01 J V* l.OOOQOD C 5 JG= 2.1000CD 00 

KP= 2. 160000 02 KI= 9.700000 03 H= 6. 000000 02 K, 0= 5 • 7 583 5D 03 

Kl= 1.37102D 03 K2= 1.000000 00 K3= 5. 000000-01 

IT ER AT I ON FREOUENC Y ~ AMPLI TUOE 

0 0.32077926 07D 01 0 . 2 57430 91 2 1 D- 05 

1 _C#96 3006 56000 00 -0.18201961810-05 

l 9. 630070-01 -1.820200-06 

I HC - 9Q0I EXECUT ION TE RM I NATING DUE TO ERROR COUNT FOR FRRQR NUMBER 26 1 


Figure 1-23c. 
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LST SYSTEM-NUMER I CRL SOLUTION OF 1 +N*GI +rUN*62= 0 . 

GRMMR= 1. 38 GOOD 07 T0= 1.000000-01 JV= 1.000000 05 JG= 2.100000 00 


KP- 2.160000 

02 K I = 9. 

700000 

03 H= 6.000000 02 

Kl= 1.371020 

03 K2= 5 . 

0 0 0 0 DO 

0 0 K3= 1 .00 0 0 OO - 03 

ITERATION 

FREQUENCY 


AML I TUBE 

0 

0 .44000000000 

01 

0 . 3 8 0 0 0 0 0 0 000- 0 5 

1 

0.43622079060 

01 

0.37386218260-05 

o 

0 . 4 3 6 6 6 2 6 9 8 6 D 

01 

0 .377 03966 090— 05 

3 

0.4371 3786550 

01 

0.37966647630-05 

4 

0.4371 3793860 

01 

0.37969710200-05 


ND= 1 FREQUENCY- 4.37138D OORAD-'SEC AMPLITUDE= 3.797390-06 NIT* 4 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N+G1 +N*N*G2= 0 . 

GAMMA* 1.380000 07 TO= 1.000000-01 JV= 1.000000 05 JG= 2.100000 00 
KP= 2.160000 02 K I = 9.7 000 00 03 H= 6.000000 02 K 0= 5.758350 03 
Kl = 1.371020 03 K2= 5.000000 00 K3= 1.000000-02 


ITERATION 

FREQUENCY 



AMPLITUDE 

0 

0.43713786550 

01 

0. 

37973907150-05 

1 

0.44394497190 

01 

0 . 

392277 164 ID- 05 

C 

0.43948472930 

01 

0 . 

38 066764 070- 05 

o 

0.43516261430 

01 

0. 

37150209610-05 

4 

0.43081 0833 OP 

01 

0 . 

-;6 0 06 188890- OS 

5 

0 .4326 02 00530 

01 

0. 

36707418870-05 

6 

0.432379439 ID 

01 

0 . 

36843250610-05 

"7 

1 

0 .4328135618 If 

01 

0. 

3 6 1 2 5 9 2 3 910- 0 5 

r» 

© 

0.43273652480 

01 

0. 

36467 079450- 05 

9 

0.43371450590 

01 

0. 

36940378420-05 

10 

0.434 06’9698i Op 

01 

0. 

37153324880-05 

11 

0 .434 06655990 

01 

0.37171261240-05 


NO- 1 FREQUENCY* 4.340390 OGRAD/SEC AMPLITUDE* 3.718220-06 NIT= 11 


Figure l-24a. 
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LST SYSTEM-NUMER ICflL SDL UT I Oh OF l+N*Gl+N^N*b£=0 . 

GAMMA® 1 .38 GOOD 07 TQ= 1.00 GOOD- 01 J'v'= 1.0 0 GOOD 05 JG= £.100000 00 
K P= £.160000 08 K I = 9. 70 GOOD 03 H= 6.0 000 OD 08 K0= 5.75335B 03 
K 1 = 1 .371 08 B 03 K£= 5.0000 OB 0 0 K3= 1 .0000 OD- 0 1 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .434 039374 IB 

01 

0.371 8883786D- 05 

1 

0 .4 19008599 ID 

01 

0 . 335339 1 99 OD- 05 

£ 

0.4 0443 1549 IB 

01 

0 . 8984578933D- 05 

3 

0 .37477 09579B 

01 

0. 83 0 0337888B- 05 

4 

0 . 38849551 59B 

01 

0 . 8566373 0 9 3 B — 05 

5 

0 . 389479365c B 

01 

0 . £5837 07884 B- 05 

6 

0 .339857c 1 8 1 B 

01 

0 .8668 049 1 68'B- U5 

7 

l 

0 . 3855 0£ 0759B 

01 

0.8538341 189D- 05 

8 

0 .37£ 158771 OB 

01 

0 « ££387643 05D— 05 

*5 

0 .38 073587 03D 

01 

0 . £39433£8£ 1 B- 05 

10 

0.388 08 035 16B 

01 

0 .£411 07884 OD- 05 

11 

0.38£fe8549£7D 

01 

0 . £4 1 88 1 33£5B— 05 

1£ 

0 .33335 051 68 B 

01 

0 .£44 153 3335 D- 05 

13 

0 .3888836738B 

01 

0 . £494 038 04£D- 05 

14 

0 .381 08 07139B 

01 

0.84837 09 089D- 05 

15 

0.447876948 IB 

01 

0 . 37454 1 6 1 46B- 05 

16 

0 .4305404333 B 

01 

0.34347071 18D-05 

17 

0.41 55898 1 88 B 

01 

0 . 3 1 34 1 7£ 1 96B- 05 

18 

0 .40 089479 1 6B 

01 

0. £7936375 19D- 05 

13 

U .663187 8 8 8 4 B 

01 

0 .943543874 IB- 05 

£0 

0.546 08 0388 ID 

01 

0.7181 84848 OB- 05 


NE= 1 FREQUENCY® 5. 46 03 OB OORAB^SEC AMPLITUDE® 7.1£lS5D-06 NIT® £0 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N+61 +N*N*b£= 0 . 

GAMMA® 1 .3S000B 07 TD= 1 . 0 000 OB- 01 JV= 1 . 00 00 OB 05 Jb= £.1000 OB 00 
KF’= £.16 00 OD 0£ KI= 3.700 0 OB 03 H= 6. 0000 OB 0£ K0= 5 . 75835 B 03 
Kl= 1.371 OED 03 K£= 5. 0000 OB 00 K3= E.OOOOOD-Ol 


ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0.546 08 0388 IB 

01 

0 .71 8 184848 OB- 05 

1 

0.4 77449746 ID 

01 

0 .49948 1 7969B- 05 

c. 

0.4351 151367B 

01 

0 .3777846 186B-05 

o 

0. 3983387 133D 

01 

0. 88855 19788B- 05 

4 

0 .86179 8 5 4 8 6 D 

01 

- 0 . 395733585 OD- 06 

£ .617 

99D 00 -3.957 

'34D- 

07 

c 

0 . 1 0 0 0 0 0 088 1 D 

01 

0 . 9 9 9999781 1 B 0 U 

s 

n . l OOOOOOOOOD 

01 

0 . 1 0 0 0 0 0 0 0 0 0 B 0 1 


HB= 1 FREQUENCY® l.OOOOOB OORHD/SEC AMPLITUDE- l.OOOOOB 00 NIT® 6 


Figure 1 -24b. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+H^1+N^G£ = 0. 

GAMMA* 1 .38 GOOD 07 TO* 1 . OOOOOD-Ql J‘v'= 1.00 0 OOP 0*.' -Jb= £.10 0 0 0 0 00 

KP= £.160000 Oc KI= 9.700000 03 H= 6.00000D 0£ KO— 5.753350 03 


* 1.371030 

03 K£= 5 . 

0 0 0 0 00 

00 K3= 1.000000 

ITERATION 

FREQUENCY 


AMPLITUDE 

0 

0 .33000000000 

01 

0 . 1 3 0 0 0 0 0 0 0 0 0 - 05 

1 

0 . 3£8 1 4b’ 1 5 9 £ 0 

01 

0 . 1344883 0380- 05 

£ 

0.33913300540 

01 

0.13539977370-05 

3 

0.33931 ££41 90 

01 

0 . 1 3536 £ 03830- 05 

4 

0. 339338479 OD 

01 

0.13543039590-05 

5 

0.33936035550 

01 

0 . 1 3543763330- 05 


ND= 1 FREQUENCY* 3. £33810 GORAD^SEC AMPLITUDE* 1 .£54980-06 NIT* 5 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N*G 1 +N+N+G3* 0 . 

GAMMA* 1.330000 07 T0= 1.000000-01 JV* 1.000000 05 JS= £.100000 00 
KP= £.160000 0£ K I = 9.700000 03 H= 6.000000 0£ K0= 5.753350 03 
Kl = 1 .371030 03 K£= 5.000000 00 K3= 1.000000-03 


ITERATION 

FREQUENY 


AMPLITUDE 

0 

0.339331 093 30 

01 

0.1 £549339470- 05 

1 

0.33017731750 

01 

0 .1368 0766 0 1 0-05 

3 

0.33064941970 

01 

0.13770117840-05 


0 .33073373510 

01 

0 • 1 £93558394 0— 05 

4 

0.33156357810 

01 

0.13005530840-05 

5 

0.33161933980 

01 

0.13013796060-05 


ND= 1 FREQUENCY* 3.316600 OORRD/SEC AMPLITUDE* 1.303330-06 NIT* 5 


Figure l-25a. 
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LST SYSTEM-NUMERICAL SOLUTION OF 1+N^G 1 +N«*N^G£=0 . 

GAMMA* 1 .38 GOOD 0? TO* 1 . 0000 OB- 0 1 JV= 1 - 000 DOD 05 JG= 2.100000 00 

KP= £.16 0 0 OB 02 K I = 9.70 0 0 OD 03 H= 6 . 0 0 0 0 OB 0£ K 0* 5 . 758350 03 

Kl = 1.371020 03 K£= 5.00 00 011 0 0 K 3= 1.000 GOD- 01 

ITERATION FREQUENCY AMPLITUDE 


0 

0 . 33 165979S2D 

01 

0 . 13 02379 026D- 05 

1 

0.34133369220 

01 

0.1 5 0448 033 ID- 05 

' L. 

0.34971371660 

01 

0 . 1 6752095380-05 

3 

0 .365 098324 7D 

01 

0.1 363445424D- 05 

4 

0 . 382258 558 1 D 

01 

0 . 233 3 5 38259 D- 05 

5 

0.4 02535 9 796 D 

01 

0.23692260760-05 

6 

0 .43 U 6719 9 6 3 D 

01 

0.360 0958328 D— 05 

7 

0 .4 163884584 D 

01 

0.3293571 0130-05 

3 

0.401 03 76 £22 D 

01 

0 . 29 038489720-05 

9 

0.3552 083 0230 

01 

0. 13745653790-05 

ib 

0 .3642164 398 D 

01 

0.20315125410-05 

u 

0.368315 13760 

01 

0.211 00277830-05 

12 

0.37050093930 

01 

0.21339590980-05 

13 

0. 37 H 63922 ID 

01 

0.209 06859450-05 

14 

0 . 379 0756437D 

01 

0.23 192 094480-05 

15 

0.381 04 0044 6 D 

01 

0 . 23356 039340- 05 

16 

0.33 14 4553320 

01 

0.25111 166800-05 

17 

0 .39939 03443D 

01 

0 . 286 62228 750- U-5 

18 

0 . 363 161 689 00 

01 

0.20519668410-05 

19 

0.37 057 053 16D 

01 

0.21761389410-05 

£0 

0.3752 088 0650 

01 

0.22536314250-05 


NO* 1 FREQUENCY* 3.752090 OORAD/SEC AMPLITUDE* 2 .253630 * 06 NIT= £0 


LST SYSTEM-NUMERICAL SOLUTION OF 1 +N<H3 1+N*N*G£=0 . 

GAMMA* 1 .33 GOOD 07 TO* 1.00 GOOD- 01 JV= 1.00 GOOD 05 Jb= 2.100000 0 0 

KP= £.1600 00 02 K I = 9.700 0 OD 03 H= 6.00 0 0 OD 0£ K 0= 5 . 75835B 03 

K 1 = 1 .371 02D 03 K2= 5 . 0 0 0 0 OD 0 0 K3= £ . 0 0 0 0 OD- 0 1 

I TERAT I DM FREQUENCY AMF’L I TUBE 

0 0.37 5 £ 0 3 8 0 6 5 D 01 0.22536314250-05 

1 0 .4209056 132D 01 0 .32729271640-05 

£ 0 . 3834 035 1 9£B 0 1 0 . £56 0546 1 35D- 05 

3 0.538 06 08335 D 01 0 . 6 1 52 1 3S662D- 05 

4 0 . 4 758 73 8 63 7D 01 0 . 46788 1 5369D- 05 

5 0 .434273 1 489D 01 0 . 368 1 723091 D- 05 

6 0 .398000471 5D 01 0 .23157169240-05 

7 0. £54453876 OD 01 -0 .56433575870-06 

1 2.544540 00 -5. 643860-07 

3 0 . 999999877 1 D 0 0 0.1 0 0 0 0 0 0 1 £50 0 1 

9 0.1 0 0 0 0 0 0 0 0 OD 01 0.10000000 0 OD 0 1 

NO* 1 FREQUENCY* 1.000000 OORAD/SEC AMPLITUDE* l.OOOOOD 00 NIT* 9 


Figure 1 -25b . 
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